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ABSTRACT 

This report presents some theoretical studies on the time -independent 
and oscillatory combustion of nonmetallized AP/composite propellants. The 
study has for its aim a coherent and unified interpretation of the voluminous 
data available from experiments related to propellant combustion. Three 
fundamental hypotheses are introduced: the extent of propellant degradation 
at the vaporization step has to be specified through a scientific criterion; the 
condensed-phase degradation reaction of ammonium perchlorate to a vapori- 
zable state is the overall rate-limiting step; gas-phase combustion rate is 
controlled by the mixing rate of fuel and oxidizer vapors. In the treatment 
of oscillatory combustion, the assumption of quasi- steady fluctuations in the 
gas phase is used to supplement these hypotheses. In comparison with experi- 
mental data, this study predicts several of the observations including a few 
that have remained inconsistent with previous theoretical results. 

As a prelude to propellant combustion studies, the behavior of AP is 
pursued in some detail. Theoretical predictions of the linear regression rates 
of AP, including explicitly the condensed-phase Arrhenius degradation term, 
are seen to match well with experimental hot plate data. Based on available 
experimental evidence of the existence of a melt layer on the surface of self- 
deflagrating AP, it is assumed that pres sure -dependent condensed-phase 
degradation in the melt layer controls the deflagration rate. The results ob- 
tained, by specifying the extent of degradation (at the vaporization step) through 
the vapor pressure equilibrium criterion, are found to predict the linear re- 
gression rate, the pressure index n, and the initial temperature sensitivity 
close to experimental data. The pressure index n is revealed as a composite 
quantity incorporating the component effects of condensed-phase degradation 
rate sensitivity to pressure and the vapor pressure effect. 

The analysis is generalized to the combustion of composite propellants, 
where the site of the rate-limiting degradation reaction of AP is assumed to be 
a thin layer on the AP particles in the propellant. The case of surface reac- 
tions in a melt layer augmenting subsurface reactions is also considered with 
a view to include propellants with readily melting binders. The temperature 
at the wall (i. e. , the interface plane between the condensed phase and the vapor 
phase) is required to be specified through proper matching with the gas-phase 
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energetics. Although a complete solution to the gas phase has remained 
elusive because of our lack of understanding of the non-laminar fluid dynamics 
encountered, two plausible models are analyzed: constant wall temperature, 

and uniform combustion in the gas phase. It is shown that either of these 
models, used in conjunction with the condensed phase analysis, yields results 
(linear regression rate versus pressure, flame standoff distance, etc. ) close 
to experimental trends, thereby de- emphasizing the importance of gas -phase 
details in the combustion of propellants. 

The response function of composite propellants, including explicitly 
the pressure-dependent degradation term in the condensed phase, is theoreti- 
cally derived. The method of inner and outer expansions with the reduced ac- 
tivation temperature parameter 0 a X= (E/R.T w ) • (1 - T q /T^) as the singular 
perturbation parameter has been applied to the problem of oscillatory com- 
bustion. It is seen that greater difficulties are associated with the specification 
of the boundary conditions than with the actual solution procedure. Two physi- 
cal situations are considered. In one case, the condensed-phase reactions are 
treated as taking place wholly in the subsurface region; in the other, surface 
reactions in a melt layer augment the subsurface degradation reactions. In 
both cases, the theoretical expression for the response function has the wall 
temperature fluctuation (f ) as an unknown quantity. Two models are con- 
sidered for the gas phase in order to determine the complex amplitude of the 
wall temperature fluctuations. In one model, the "flame" temperature fluctua- 
tions are related to chamber pressure fluctuations through the isentropic rela- 
tion. (It is shown that this assumption, however, does not lead to isentropic 
fluctuations at the wall. ) In the other model, the gas phase combustion proc- 
esses are assumed to be uniform even during oscillatory combustion. The re- 
sponse functions, so derived, exhibit dependence on mean chamber pressure; 
strongly so with the model of adiabatic fluctuations in the gas phase and weakly 
for the case of uniform combustion. For the models of no melt layer and non- 
oscillating melt layers, the theoretical response functions are seen to ap- 
proach large values at low frequencies (as observed in some experiments) with 
the adiabatic assumption for the gas phase fluctuations. The proper limit (i. e. , 
the steady-state pressure index n) is, however, reached at zero frequency 
when it is recognized that very slow changes in the gas phase are isothermal 
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and not adiabatic. It is also seen that a "zero-n" propellant can exhibit 
fairly strong instability behavior. 

Thus, it is seen that many of the apparently diverse experimental ob- 
servations are all consistent within the framework of the present theoretical 
developments. It is, therefore, concluded that the results obtained so far 
are sufficiently encouraging to warrant further research on similar lines. 
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NOMENCLATURE 

constant in the uniform combustion law, equation (72) 

r -3 -l -l-i i 

Lgm. cm. sec. atm. J 

mean size of the oxidizer particles in the composite propellant 
[cm. ] 

pre- exponential factor in the Arrhenius law for thermal degrada- 
tion [sec. ^ 1 \ 

constants (with i = 1, 2, 3, ... ) used as shorthand notations to 
simplify the algebra [dimensionless "j 

constants (with i = 1, 2, 3, ... ) used as shorthand notations to sim- 
plify the algebra in the treatment of liquid layers on composite 
propellants [dimensionless] 

specific heat [cal. gm. " * °C ^ 1 

heat of degradation of the solid, i. e. , heat required to convert 
one gram of the polymer or crystal into one gram of the individual 
repeating units [cal. gm. ~ * ] 

mutual diffusion coefficient for the oxidizer and fuel gases 
Lem. sec. j ! 

activation energy for thermal degradation [cal. mole” * ] 

statistical mean fragment size vaporizing [dimensionless] 

statistical mean fragment size at the solid-liquid interface 
[dimensionless] 

normalized heat of degradation [dimensionless] 

coefficient of thermal conductivity [cal. cm.’* sec.’* °C’*] 

thickness of the surface melt layer in composite propellant 
combustion [cm. ] j 

thickness of the surface melt layer in the self-deflagration of AP 
single crystals [cm. ] 

molecular weight [gm* (gm. mole )” * ] 
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NOMENCLATURE (cont'd. ) 

r -2 - 1 , 

mass flux of combustion gases Igm. cm. sec. J 
mole fraction of species i [dimensionless”] 

gas phase combustion rate during uniform combustion 

r " 3 -In 

Lgm. cm. sec. ] 

fractional number of backbone bonds referred to the number in 
the undegraded state [dimensionless"] 

empirical index of pressure in the steady burning rate law for 
propellants [dimensionless] 

pressure [atm. ] 

normalized temperature gradient [dimensionless] 

heat released by combustion gases upon complete combustion 
[cal. gm." 1 ] 

Reynolds number for the gas phase processes [dimensionless 1 

constants (with i = 1,2,3,...) used for shorthand notations 
[dimensionless] 

universal gas constant [cal. mole * °C *] 

complex (pressure) response function [dimensionless] 

linear regression rate [cm. sec. ] 

normalized heat release rate in the gas phase (defined in equation 73) 

[ dimens ionle s s ] 

temperature [°K] 

time coordinate [sec. 1 

shorthand notations (with i = 1, 2, 3, . . . ) us ed in the section on 
uniform combustion model for the gas phase [dimensionless] 

gas (mass flow) velocity above the burning propellant [cm. sec. - *] 

volume fraction of oxidizer in a composite propellant [dimensionless] 
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NOMENCLATURE (cont'd. ) 

normalized temperature difference between the flame zone and 
the wall [dimensionless 1 

distance coordinate [cm. "| 

normalized distance coordinate [dimensionless] 

shorthand notations (with i - 1, 2) used in the treatment of 
oscillatory combustion [dimensionles s 1 

empirical constant in the FSV rule [atm.^] 

empirical index of pressure in the FSV rule L dimensionless ] 
ratio of specific heats of combustion gases [dimensionless] 
empirical constant in the FSV equation [°K * ] 

a small parameter used in the expansions; any first order quantity; 
r'/r [dimensionless] 

normalized flame standoff distance [dimensionless] 

normalized stretched temperature coordinate in the "inner" 

region in the asymptotic analysis of degradation reactions [dimensionless 

activation energy parameter, = E/BT w [dimensionless] 

thermal diffusivity [cm. sec.” ] 

normalized mass burning rate in the gas phase [dimensionless] 

normalized time-independent regression rate eigenvalue 
[dimensionless] 

normalized time -dependent regression rate eigenvalue 
[dimensionless] 

complex root of the heat transfer equation [dimensionless] 

normalized amplitude of pressure fluctuations [dimensionless] 

normalized "inner" dependent variable in the asymptotic analysis 
of degradation reactions [dimensionless] 
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NOMENCLATURE (cont'd. ) 

_ 3 

density [gm. cm. ] 

normalized temperature [dimensionless] 

thickness of the surface melt layer on the oxidizer crystals in a 
composite propellant [cm. 1 

\ 

normalized temperature parameter used in solid phase analysis, 

= (T -T )/ T [dimensionless] 
wow 

—2 

normalized frequency of fluctuations, = urn/ r [dimensionless] 
frequency of fluctuations [sec."'] 
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plane of burning, i. e. , flarnie zone 

freezing (of chemical reactions) plane in the solid; flame region 
in the gas 

gas 

deep solid (ambient) conditions 
solid 

solid-liquid interface plane 
subsurface region 
wall plane 

outer region in the asymptotic analysis 

time- averaged part 

fluctuating part 

time -dependent quantity 

reference quantity 



SOME CONSIDERATIONS IN THE COMBUSTION OF 


AP/ COMPOSITE PROPELLANTS 


I. INTRODUCTION 

Solid propellant rockets are here to stay. It is the wish of re- 
searchers in the field to be able to say the same of propellant combustion 
theories. Continuing development of solid propellant rockets has there- 
fore depended largely on experimental programs which have proved ex- 
tremely useful so far as the technological application is concerned. How- 
ever, the lack of quantitative correlations among the vast amount of ex- 
perimental data has left the predictive design of rockets an unrealized 
ideal. Many deficiencies in our understanding that are adequately con- 
cealed by empiricism when one considers steady combustion are revealed badly 
when the propellant enters unsteady and oscillatory modes of combustion. 

Since it has not been possible at the present time to accurately predict 
whether a given propellant would burn in the steady mode or in the unsteady 
mode in the rocket chamber, the difficulties are easy to comprehend. 

Even a coherent and unified interpretation of experimental data has re- 
mained a challenging problem, so that the few isolated successes of theories 
have not been above scepticism. Thus, the need is felt for work at the fun- 
damental level in an attempt to identify the core processes common to 
most, if not all, of the propellants in use. If successful at the primary 
task of coherently interpreting different experimental phenomena in terms 
of a few unifying concepts, the study could then work out the details and re- 
finements necessary in the prediction of results of any specific experiment. 
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In the present work a study is made of ammonium perchlorate - 
based, non -met sill ized composite propellants. The objective is to identify 
the fundamental processes that are likely to be common in a variety of pro- 
pellants and propellant applications. The motivation for the present work 
comes from the observation that under normal conditions the rate processes 
in the gas phase are likely to be much faster than those in the condensed 
phase. As elaborated on in Section II, the present work differs from those 
available in the literature in three important aspects. The chemical kinetic 
degradation reactions in the condensed phase are explicitly included in the 
analytical treatment. The importance of a scientific criterion in specifying 
the extent of propellant degradation before vaporization is stressed, and it is 
shown that the usual arbitrariness is removed by applying the vapor pressure 
equilibrium criterion at the propellant surface. Lastly, the gas phase chemi- 
cal reactions are treated as wholly controlled by pressure-independent mo- 
lecular mixing processes. It is found that the study predicts several of the 
features observed in AP/propellant combustion. 

Although no fundamentally new concept is introduced, some of the 
thoughts are either relatively new to, or have not all been considered at the 
same time in, propellant combustion theories. Those aspects of the problem 
that have already received adequate treatment on similar lines are not pur- 
sued here. Also, no attempt is made to review the literature on propellant 
burning. Not only is the research field very active at the present time, but 
the literature on the subject is also very extensive, rendering it difficult to 
cite a few references in the limited space here. Nevertheless, the reader 

should find useful information on the general subject in the collection by 

27 

Warren , on the work of the Princeton group in Steinz, Stang and Summer- 

11 v 28 

field , on the detailed theoretical studies in Williams, Barr'ere, and Huang , 
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and of the Russian work in Bakhman and Belyaev and Novikov, Pokhil, and 
30 31 

Ryazantsev . The monograph by Price and Culick familiarizes the reader 

with the status of the problem in 1969, which has not changed much since then. 

After the present study was initiated, it was found that some of the thoughts 

run parallel to those of Willfred Schmidt, who has also been considering the 

19-21 

importance of condensed phase reactions in his recent studies 

Since AP is a prominent constituent, its behavior by itself is studied 
in some detail in Section III. The condensed-phase degradation rate obtained 
from small- sample isothermal data and the fragment size of the AP parti- 
cles leaving the wall obtained from a vapor pressure criterion are used to 
predict the overall linear regression rates of AP in different experiments. 
Very good agreement is found with hot-plate data and with single- crystal 
deflagration experiments. It is found that the pressure index n is actually 
a composite quantity incorporating both the effects of pressure-dependent 
chemical reactions and the vapor pressure effect on the fragment size va- 
porizing. It is also found that a power law for pressure dependence of re- 
gression rate is an approximate representation at best. The actual depend- 
ence is non- simple and involves at least one logarithmic variation. The 
effects of initial temperature variation on the linear regression rate are also 
predicted and are found to be close to experimental observations. 

Time -independent burning of AP composites is studied in Section IV. 
The study is undertaken in the belief that a thorough understanding of the os- 
cillatory combustion of composite propellants would be difficult in the absence 
of an under standing of the time -independent burning. Effects of burning rate 
catalysts and oxidizer particle size are explicitly included. The depolymer- 
ization reactions of fuel (binder) are also included. The general problem of 
gas phase details is formulated and a solution is written in terms of non- 
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dimensional groups, which are also recognized to be important similarity- 
parameters. A completely self-contained solution to the gas phase has not 
been possible because of our lack of understanding of the gas phase fluid 
dynamics above the regressing surface of the propellant. In order to render 
the solution self-contained, two physical models are considered; constant 
wall temperature and uniform combustion in the gas phase. It is seen that 
either of the models can predict results close to experimental trends (both 
qualitatively and quantitatively). This reinforces the belief that precise de- 
tails of the gas phase mechanics are not very important in determining the 
general trends, if the rate-limiting reactions occur in the condensed phase. 

Oscillatory combustion of composite propellants is considered in 
Section V. Employing the method of "inner" and "outer" expansions, linearized 
analytical solutions are obtained for the response function, including explicitly 
the Arrhenius reaction rate term in the condensed phase. Two physical situa- 
tions are considered. In one case, the condensed phase reactions are treated 
as taking place wholly in the subsurface region; in the other, surface reactions 
in a melt layer augment the subsurface degradation reactions. In both cases, 
the theoretical expression for the response function has the complex amplitude 
of wall temperature fluctuation as an unknown quantity. Two models are con- 
sidered for the gas phase (treated quasi- statically ) in order to determine the 
complex amplitude of the wall temperature fluctuations. In one model, the 
"flame" temperature fluctuations are related to chamber pressure fluctuations 
through the isentropic relation. (It is shown that this assumption, however, 
does not lead to isentropic fluctuations at the wall. ) In the other model, the 
gas phase combustion processes are assumed to be uniform even during oscil- 
latory combustion. The response functions exhibit dependence on mean 
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chamber pressure; strongly so with the model of adiabatic fluctuations in the 
gas phase and weakly for the case of uniform combustion. It is also seen 
that a "zero-n" propellant can exhibit fairly strong instability behavior. 

The present study is not complete. However, the results obtained so 
far are thought to be encouraging. The basic conclusion at this stage is that 
it impossible to understand a variety of superficially different experimental 
phenomena within the broad framework of a few unifying concepts. 
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II. QUALITATIVE CONSIDERATIONS 

2. 1 Introduction 

The fundamental physical structure of the present work is discussed 
in this section. It is also the aim here to anticipate, on physical grounds, 
the analytical results of subsequent sections. The models of condensed 
phase reactions, the concept of fragment size vaporizing (FSV), the gas 
phase processes, and a brief discussion on the numerical values of the 
thermophysical and thermochemical constants that are needed later are in- 
troduced separately. 

2.2 Condensed Phase Details 

In an attempt to trace the entire history of the oxidizer (or fuel) 
from the deep solid state to the product state, the condensed phase details 
arise in a natural way. Since many composite propellants are heavily loaded 
with ammonium perchlorate (80 per cent or higher), we study the AP be- 
havior at first. In its unaffected state AP may be looked upon as a large 

"molecule" made up of the fundamental building blocks -(NH^C-tO^)- . 

7 12 13 

While the degradation of AP has been studied in detail ’ ’ and is known 

to involve the production of ammonia (NH^) and perchloric acid (HC-fD^), it 
is interesting to inquire into AP behavior prior to such a decomposition. It 
would seem logical to expect that in the sequence of decomposition reactions 
the first stage would be the scission of the weak bonds between the neighbor- 
ing -(NH^C-tO^)- units, particularly when the rate of heating is very high, 
as in propellant applications. The large crystals degrade into smaller 
groups of -(NH^C-tO^)- and we ask ourselves whether the degradation of AP 
completely into single molecules of (NH^Cf,0^) and the decomposition to 
NH^ and HC70^ are prerequisites for the oxidizer to leave the condensed 
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phase and enter the vapor phase. It would seem unlikely that such indeed 

is the case. Under certain conditions, pure AP has been known to sublime 

14 

and leave the surface as pure AP. In one extreme case, at very low (sub- 

15 

atmospheric) pressures, macroscopic particles of AP have been observed 
to get ejected from the surface (of a burning composite propellant) into the 
vapor phase. These observations suggest that in order to specify, as a 
function of chamber pressure and wall temperature, the actual state of the 
oxidizer leaving the surface, a scientific criterion is necessary. The prob- 
lem then is almost identical with that associated with the combustion of any 
solid, like a polymer, for example. Specification of the size of AP at the 
surface requiring that the vapor pressure sum of all fragments leaving the 
surface at the wall temperature equal the chamber pressure appears to be 
a scientific criterion. This vapor pressure criterion is invoked in the 
present study. 

The fragment size of AP leaving the wall, as a multiple of the funda- 
mental unit - (NH^C-fO^)- will be designated the fragment size vaporizing, 
or FSV, for short. For a constant wall temperature, the FSV will be very 
large at low pressures and small at high pressures. At sufficiently high 
pressures, the fragment size specified by the vapor pressure criterion gets 
smaller than the fundamental unit (NH^C^O^) . Decomposition into smaller 
molecules is anticipated in such cases. In any attempt to experimentally 
determine the actual fragment size at the wall, it should be remembered 
that the entire process of propellant regression is a nonequilibrium one and 
that rapid quenching of all reactions is necessary, immediately after the 
species enter the vapor phase, in order to study the species leaving the 
surface. That is, reactions (degradation, decomposition) in the vapor phase 
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could mask the identity of the species actually leaving the wall surface. 

At the surface of a burning composite propellant, both fuel and oxi- 
dizer species are present. Thus, the necessity for a proper mixing rule 
arises in the generalization of the vapor pressure criterion to multicompo- 
nent equilibria. A detailed study would consider the energetics as well and 
include the relative strengths of bonds between the neighboring molecules 
in the fuel and oxidizer in the propellant, the mole fractions of fuel and 

li 

oxidizer and separate wall temperatures, if the "two-temperature" concept 
is a physical reality. The nature of the present work does not require a 
sophisticated treatment of these effects. A logical mixing rule in Section 
IV handles the problem adequately. For a thorough discussion of the vapor 
pressure criterion, the reader is referred to ref. 17. 

In order to determine the numerical value of the FSV as a function 
of specified conditions, a study has been made of the vapor pressure data 
of hydrocarbons and a rule has been evolved in Appendix B to predict the 
FSV. The approximate nature of extrapolations to AP and polymer from 
hydrocarbon data is recognized, but it is felt that for a simple physical 
quantity like the vapor pressure, the chemical nature of the molecules is 
not crucially important. As a matter of fact, the experimental vapor pres- 
sures of hexane and methylmethacrylate monomer (identical molecular 
weights, 100) are very close over a range of temperature*. 

A non-integer value of FSV predicted by equation B-2 is to be under- 
stood as an average over all fragment sizes and not as the presence of 
breakdown products from the fundamental repeating units. 



-9- 


From the above discussion, there emerge the following three con- 
siderations which form the backbone of the rest of this study. 

(i) The fragment size of vaporizing AP can be specified through 
vapor pressure equilibrium criterion. 

(ii) The degradation of macroscopic crystals of AP into vaporizable 
fragments of AP is the fundamental rate-limiting reaction of 
interest. 

(iii) Assuming that the AP degradation is brought about by reactive 
species, the degradation rate is taken as directly proportional 
to pressure. 

It is seen that the above three considerations not only remove con- 
siderable arbitrariness in propellant combustion theory but also lead to re- 
sults that are close to experimental trends both qualitatively and quanti- 
tatively. 

2. 3 Gas Phase Details , 

It is useful to have an order -of-magnitude estimate of the character- 
istic scales in the physical problem. Chemical kinetic rates in the gas 
phase may be inferred through measurements of aerodynamic flame speeds 
in experiments where the speed is known to be controlled by the chemical 
kinetic reaction rates, such as premixed laminar flames. At atmospheric 
pressures, typically the aerodynamic flame speeds of air /hydrocarbon 
mixtures are like 100 cm/sec. The characteristic length scale, which is 

the flame thickness, is of the order of 100 |_i. Hence, the characteristic 

-4 

chemical reaction time is like 10 sec. At high pressures, as in propel- 
lant applications, we expect the reaction time to be much smaller. If the 

_ 5 

reaction is bimolecular, the characteristic time at 10 atm. would be 10 sec 



- 10 - 


The characteristic transport/mixing time in the burning of a com- 
posite propellant is at most of the order of the ratio of combustion zone 

standoff distance to the mean velocity of the gases leaving the propellant 

-4 

surface. Employing typical numbers we arrive at a time like 10 sec. for 
the mixing process. While these estimates are not conclusive, the general 
trend is in favor of transport/mixing control and not of chemical kinetic 
control. Another factor that greatly accelerates the chemical kinetic rates 
in propellants is that the oxidizers used are much more powerful than air. 
Perhaps the most direct support to the argument is provided by the obser- 
vation that the actual high reaction-rate zone occupies a small fraction of 
the totcil "flame" standoff distance at conventional pressures. 

Thus, we look for pressure-dependent mechanisms other than gas- 
phase chemical kinetics, leading to the examination of condensed phase de- 
tails mentioned in II. 2 . However, we do need a solution to the gas phase to 
completely specify the regression rate. 

The gas phase details have received a great deal of attention over 
the years. Several models have been presented. Actually, when once the 
assumption is made that the gas phase is not rate controlling, any model 
that incorporates the essential heat transfer details will suffice. This is 
because the role of the gas phase above the burning propellant degenerates 
from an active control of general propellant behavior to one of supplying 
boundary conditions on the solid phase details. 

The reactants burn in a non-premixed combustion zone in the gas 
phase except in the following two cases. 

(i) Gas phase chemical -kinetic rates become very slow, because 
of very low pressures or special ingredients, so that molecular 
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mixing processes take the gases to a premixed state before 
combustion. 

(ii) A thorough mixing of fuel and oxidizer takes place in a surface 

1 8 

layer on the propellant before they enter the vapor phase. 

The details of the gas phase, as required in the solid phase analysis, 
are developed in Section IV. It is seen that the simplest model of the com- 
bustion zone (treated like a "black box") is adequate for our purposes. The 
question of laminar versus turbulent nature of the fluid process is discussed. 
The important parameters are identified. A formal solution is obtained, 
although a closure has not been possible, mainly because of our lack of 
understanding of transport and mixing processes in propellant burning. 

2. 4 Numerical Values of Constants 

Although the number of fundamental quantities that are needed in 
the present study is not large, considerable uncertainty exists in the nu- 
merical values of even the few properties that find extensive use. The prin- 
cipal quantities are the thermal diffusivity of the propellant, the pre- 

/ 

exponential factor and the activation energy of the degradation reaction 
and the heat of degradation. Early in the present study it was decided to 
use a single set of values consistently through the work. 

(i) Thermal Diffusivity, x . A majority of polymers have a thermal 
- 3 2 

diffusivity around 10 cm /sec. , although there are indications that the 
value decreases considerably at high temperatures. AP is crystalline, 
and one would expect its thermal conductivity to be much higher than that 
for an amorphous polymer. For a material like composite propellant grain, 
we would expect the thermal diffusivity to be between that for the base fuel 
and that for crystalline AP. The values used here are: 




-12 


pure polymer base 

10“ 3 

normal AP composites 
with heavy AP loading 

1. 1 x 10" 3 

pure AP 

1. 5 X 10" 3 


Whenever a specific propellant is considered, if the h for that propellant 
is available, such a value is used in preference to the above table. 

(ii) Fundamental Rate Data for Pure AP. The high temperature 
(isothermal) degradation values of ref. 7 (page 41) are used. The values 
of the smadlest sample are used for obvious reasons. 

B = 9. 2 x 10 7 sec -1 , 

E = 28. 9 k cal/mole. 

(iii) Heat of Degradation, D. The heat of degradation is taken as ap- 

1 3 

proximately 585 cal/gm for AP. The same value is used for AP-based 
composites also. It is also taken as endothermic. While the overall py- 
rolysis of AP may be exothermic, it is possible that the rate-limiting degra- 
dation reaction is endothermic. A simple bond -breaking degradation scheme 
would suggest endothermic ity. Besides, the value of the linear regression 
rate of propellants is theoretically seen to be not a strong function of the 
heat of degradation (see later). Physically, this is because the bulk thermal 
sink contribution, namely c(T w -T q ) , far overwhelms the degradation heat 
term under most conditions. 
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III. DEGRADATION AND DEFLAGRATION OF 
AMMONIUM PERCHLORATE 

3. 1 Postulated Model 

The aim of the present section is the theoretical prediction of the 
linear regression rate of AP as a function of chamber pressure (or other 
experimentally determined conditions). We examine the model (Section 3. 2 
within the framework of a one -dimensional picture. 

The overall model is depicted in fig. 1. We consider the steady- 
state degradation of a semi-infinite mass of AP (all transverse gradients 
zero). In the coordinate system used, the wall surface is held stationary. 

A plane parallel to the wall surface moves up with time and we follow the 
changes in this plane as it moves from x = +co (deep solid) to x = 0 (wall). 

Its temperature increases from the deep solid value (T q ) to the wall tem- 
perature (T w ). The mean fragment size of the AP particles changes from 
a very large value ( -*oo) at x = oo to FSV at the wall (x = 0). 

Following the numerous studies in the field, the degradation of AP is 
modeled as a first-order Arrhenius reaction. The pre- exponential factor is 
taken as directly proportional to the pressure in the molten state (see later) 
and independent of pressure in the solid. 

In the analysis to follow, energy balance due to conduction, "convec- 
tion, " and heat of degradation are considered. Diffusion of small fragments 
through larger ones is neglected. That is, non-random velocities of indi- 
vidual fragments differing appreciably from the mean are not considered. 

I 

For reasons that will be clear later, the cases of subsurface and surface 
reactions are treated spear ately. 
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3. 2 Analysis of Subsurface Reactions 

Although we are interested at present in the steady state only, the 
full (time -dependent) equations are written down since these are needed 
later (Section V) in the analysis of oscillatory burning. 

Governing Equations: 

Energy: 

2 

k-^-T + cp r - pc -|y = DpNB • exp(-E/RT) (1) 

3x 

Bond Conservation: 

- ^ = NB • exp(-E/ftT) (2) 

Boundary Conditions: 

x = 0: T=T ;x=co,T=T (3) 

w o 

x = 0 : N = N = 1 - ; x = oo , N = 1 (4) 

w FSV 

In the above equations, k stands for the coefficient of thermal con- 
ductivity (cal/cm. sec°C) , T for absolute temperature (°K) , x for the 

coordinate into the solid (cm), c for the specific heat of the solid (cal/gm. °C), 

3 

p for the density of the solid (gm/cm ) , t for the time coordinate (sec), 

D for the heat of degradation of the solid, i. e. , heat required to convert 
one gram of crystalline AP into one gram of individual ( NH^C-f.O^ ) mole- 
cules (cal/gm), N for the fractional number of backbone bonds referred 
to the number at infinity, B for the pre- exponential factor (sec - *), E 
for the activation energy of degradation (cal/mole), and ft for the universal 
gas constant (cal/mole°C). When the initial number of backbone bonds is 
very large, it is easy to see that the fractional number of bonds at the wall 
is essentially (1 - 1/FSV), an approximation that has been used here. (See 
Appendix B for details of FSV. ) 
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Defining 


and 


cprx/k = rx/ ; 


t = (T-T )/ (T -T ) 
o wo 


and substituting them into equations (1) and (2), we get (now we consider 
only the time-independent case): 


d t + dT _|< DpNB exp(-E/ RT ) 


dy 


dy 


( P cr) (T w- T o } 


and 


dN , dr kNpcB • exp(-E/RT) 

dT dy . ,2 

7 (per) 

The nondimensionalizations, 

p = dT / dy 

h s D/c(T -T ) 
w o 


and 


A = kB/pcr 

transform equations (5), (6), (3), and (4) to 
pp' +p = AhNexp(-E/RT) 
pN' = AN exp(-E/ RT) 

with 

T = 0:p=0, N = 1 

and 


T = 1 : p = p 


w 


N = N 


w 


(5) 


( 6 ) 


(7) 

( 8 ) 


where a prime denotes d/dT . 

The Value of p__ . At the wall interface, the balance of ener 

' V/ 


gy 


requires 
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heat supplied 


heat used in raising 


heat used in back- 

into the 

= 

the temperature 

+ 

bone bond 

material 


from T to T 


breaking 

- 


u o w 


-J 


“ k ^ = r P c ( T w- T o ) + Dpr/FS V 


Within the framework of the model, heats needed for phase changes need 


not be included explicitly. 

In terms of nondimens ional variables defined earlier, heat balance 
at the interface takes the form 


w 


= -1 - 


FSV 


(9) 


The Value of N . We have already seen that (Appendix B) 


N w 1 " FSV 


( 10 ) 


Solutions 

The close similarity between equations (7) and (8) may be exploited 
in a manner analogous to treatments of laminar flame propagation in pre- 
mixed gases. Multiplying equation (8) by h and subtracting the resulting 
equation from equation (7), we get an equation that can be integrated in 
closed form, leading (after the use of boundary conditions) to 

N = (p + T + h)/h . ' (11) 

Substituting equation (11) into equation (7), we get 

pp 1 +p = (h+p+T) A exp(-E/RT) (12) 

T = 0 : p = 0 ; T = 1 : p = p w (13) 

Thus, we have succeeded in reducing the initial system of third- 
order differential equations to a single first-order differential equation. 
The nonlinearity is retained. 
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Typically, the activation energy E has a value like 30 k cal/ mole, 
while the maximum temperature (which is at the wall ) is like 900 - 1000°K. 
Hence, E/ftT is at least 10-15. This implies that the R. H. S. of equation 
(12) falls to exponentially small values even at short distances from the 
wall. 


Such a behavior is particularly suited for a matched asymptotic 
analysis. Essentially, we neglect the reaction term from equation (12) 
far from the wall and obtain the solution to the resulting linear equation. 
Next, we consider the region close to the wall where the reaction term is 
explicitly included but the temperature range of interest is only a first 
order quantity. That is, expansions in powers of a small parameter are 
possible in the familiar fashion. We match the solution near the wall 
region with that away from the wall and in the process determine the re- 
gression rate, r. All of the details are available in Appendix A (and can 
also be inferred from the very similar treatment of oscillatory case pre- 
sented here in Section V). The end result is 
A exp( - 8 ) 


X = 


6 a* 


(h+1) ^' n { lUrr} 


FSV 


(14) 


where 0 s E/ftT and X = (T -T )/ T 
a w wow 

In terms of dimensional variables. 


r hB exp(-9 a ) 2 

= l 9 xX ^ J 


(15) 


Equation (15) represents the desired result which will find extensive ap- 
plications hereafter. 
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3. 3 Comparison with Experimental Results 

Recalling that the derivation in Section 3. 2 did not make provision 
for surface reactions, we should examine data from those experiments 
where there is little possibility of such surface reactions. The familiar 
hot-plate experiments seem to come under this category. However, there 
are indications that a criterion other than the vapor pressure criterion is 
needed. There is no clearly-defined equilibrium interface between a con- 
densed phase and a vapor phase. It was seen * that a constant molecular 
weight of 900 appears to match experimental data very well for the polymer 
PMMA. It is possible that a mechanical strength criterion is more applica- 
ble for hot-plate experiments. That is, as the material degrades from 
within the deep solid, a plane is reached where the increasing temperature 
and decreasing physical strength force the material out of the surface. A 
molecular weight of 900 corresponds to an FSV like 7-8 for AP (funda- 
mental molecular weight 117. 5 gm/gm. mole), and a value of 8 is used in 
the present study. At these high values of FSV, 7 or 8 will not make more 
than about 7 per'cent difference in the final regression rate, an error that 
is much smaller than the general levels of uncertainty in such experiments. 

Presented in fig. 2 are the experimental data points collected by 
9 

Powling from many different sources. The theoretical prediction of 
equation (15) is also plotted. Good agreement is evident. The variation of 

regression rate with surface temperature is extremely well predicted. It 

I 

is noted that because of the square root factor in equation (15), an overall 
activation energy of nearly 15 k cal/mole would be inferred by forcing through 
the hot-plate data an Arrhenius expression. Actually, such curve fit 
procedures grossly average the fundamental processes. The valid pro- 
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cedure is to use the isothermal Arrhenius parameters (obtained from ex- 
periments on very small samples of AP ) in equation (15), as has been done 
here. As regards the quantitative agreement of equation (15) with experi- 
mental data, it is worth remembering that no more than fundamental thermo- 
chemical property values were needed for the prediction. 

While the validity of the model is borne out by the hot-plate data, 

several inconsistencies arise when an attempt is made to generalize the 

‘ 

above model to other cases -- notably self- deflagration flame-heated AP. 

The quantitative, and even qualitative, disagreements noted below lead to 
the developments in Section 3. 4. 

To examine the simplest case, let us consider first pure crystals 
of AP. There are clear evidences in the literature that the surface tem- 
perature is no more than 900 - 950°K. In fig . 2 we note that at such wall 
temperatures, a regression rate around 0. 1 cm/sec is expected. (Actually 
r would be less theoretically because the FSV at the high pressures en- 

i 

countered is much smaller than 8, a number that has been used in the pre- 
diction presented in fig. 2. ) Typical values of r observed in AP self- 

4 

deflagration flame-heated samples are like 0. 5 to 1 cm/ sec. Thus, sub- 
surface thermal degradation rate alone cannot account for the high regres- 
sion rates observed in AP self-deflagration flame-heated samples. Secondly, 
the wall temperature in pure AP deflagration is known not to vary too much 
with regression rate. While this point has been the subject of much discus- 

i 

sion.in the past, perhaps the most carefully performed experiments that 

5 

were recently made available by Maltzev provide a very direct support. 

In any case, the near constancy of surface temperature over wide ranges 
of regression rate suggests that some mechanism that does not depend on 
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the temperature effect in the Arrhenius reaction term is in operation. It 
is readily supposed that the pressure-dependent pre- exponential factor is 
involved. 

However, equation (15) predicts a pressure index n of less than 0. 5 
always (this point will be elaborated on shortly); observed values of n are 
like 0. 7 to 0. 8 . While it may be possible to account for all of the incon- 
sistencies retaining the subsurface reaction model (thus retaining equation 
(15)), it would require fairly complicated interrelated variations in the 

simple parameters. It is shown in Section 3. 4 that a very much simpler 
■ ■ . i 

picture is also capable of accounting for all of the observations. 

i 

3.4 Self- Deflagration Flame -Heated Single Crystals of AP 

Let us make the simple postulate th at surface degradation by re - 
active species aids subsurface thermal degradation in producing vaporiza - 
ble fragments and that at high pressures completely overwhelms the sub- 
surface contribution . The direct dependence of the pre- exponential factor 
on pressure, that was postulated earlier, may lead us, at first sight, to 
predict the pressure index n to be unity (since the wall temperature is 
known to be reasonably constant). However, the fragment size (FSV) goes 
down as the pressure increases at constant wall temperature, and the addi- 
tional requirement on kinetic degradation rate at the surface keeps the pres- 
sure index below unity. We shall make a numerical calculation shortly, 
but before undertaking such a calculation, it is necessary to consider the 
physics of heterogeneous reactions in some detail. 

We expect the heterogeneous catalytic degradation by reactive spe- 
cies to take place in a very thin (but finite thickness) layer on the surface. 







- 21 - 


It would seem physical to think of this surface degradation as taking place 

in a layer where the catalytic species and the AP mixed thoroughly on the 

microscopic scale. The observation of a "liquid" layer on the surface of 
4 

deflagrating AP lends credence to such a picture. A question may arise 
with regard to the rate limiting process. It is easily seen that the diffu- 
sion of the catalytic reactive species into the melt layer is a much faster 
process than the degradation reaction. Eyring notes (see ref. 6) that dif- 
fusion in liquids may be looked upon as a "chemical reaction" in which the 
weak bonds between neighboring molecules are broken by the diffusing 
species' energy. Such bonds have strengths that are obviously equal to the 

heat of vaporization of the liquid. The heat of vaporization of many such 

! 

liquids is of the order of 5 kcal/mole. Usual activation energies of degra- 
dation are around 30 kcal/mole. Naturally, diffusion is much faster than 
degradation. 

7 

Thickness of such a liquid layer on a regressing surface is observed 

g 

to, and can be theoretic ally shown to, decrease with increasing regression 

5|C 

rates at a constant wall temperature. Finally, at a high enough regression 
rate, the melt layer may occupy only a very small portion of the surface. 

(In the ideal, one -dimensional case, the liquid layer disappears completely. ) 
When this happens, the degradation by the reactive species becomes very 
slow because diffusion through the solid is a much slower process than 
through a liquid. The process of degradation is now more likely to be 
completely controlled by subsurface thermal degradation. We thus expect 
a marked drop in the regression rate after such a point is reached. An ex- 
ploratory calculation has indicated that the sudden dip observed in the AP 

regression rate (around 2000 psia) may be due to such a phenomenon. 

Also explained in Appendix C. 
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Now we are in a position to write the relevant equations. 

The linear regression rate may be written 

/degradation rate in the \_jZthe number of bonds to'j 
\melt layer of thickness IJ ’ Vbe broken 


/ 


r = 


l N B- exp(-E/RT ) 
r w 


( 16 ) 


i l " FSSlI “ i 1 “ FSv) 

Use has been made of the fact that, at any plane, the fractional number of 
bonds left is 

1 


1 


fragment size at that plane 
and FSSL stands for the mean fragment size at the interface between the 
solid and the liquid layers. 

Consistent with the general kinetics schemes of such reactions, B 
is taken as B q P, where P is the pressure and B q is the pressure-independ- 
ent reference value of the pre- exponential factor. 

7 

Below 2000 psi, surface photographs of deflagrating AP crystals 
suggest considerable activity in the melt layer. Bubbles from below seem 
to stir the layer. Consistent with this picture, we assume a uniformly 
mixed layer. The average number of bonds in the layer follows immedi- 
ately as 


1 


N = 


FSV 


+ 1 - 


1 


FSSL 


N = 1 


‘ (fsv + fssl)/' 


(17) 


FSV FSSL 

When the surface contribution to degradation rate in the melt layer far 
overwhelms the subsurface thermal contribution, FSSL is very large com- 
pared to FSV. Thus, equations (16) and (17) may be simplified to read 


r = J,N* FSV- B- exp(-E/RT ) 

r w 


( 18 ) 
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and 

N = 1 - \ FSV . (19) 

Combining equations (18) and (19) we get 

r = £(FSV-i)B o (P/P o )exp(-E/RT w ) . (20) 

Predictions of r through equation (20) are compared in fig. 3a for 
three assumed values of the melt layer thickness £ . Experimental data 
trend is also indicated. 

Again we note the extremely reasonable prediction through very 
simple modeling. A constant value of melt layer thickness (between 5p 
and 8 p ) is seen to match experimental data quantitatively. However, since 
the melt-layer thickness l decreases with increasing regression rate (in- 
creasing pressure) as shown in Appendix C, this agreement should be due to 

self-compensating effects of slight increase in wall temperature with pres- 

5 

sure (see Maltzev ) and the decrease in the melt layer thickness. 

3. 5 Effects of Initial Temperature (T ) yariation 

^ I 

Equation (15) predicts the steady- state regression rate as a function 
of the temperature at any plane and the mean size of the fragments at the 

I 

same plane. Thus, we see through fig. lb that we can predict the regression 
rate of single crystals of self -deflag ration -flame heated AP if we use T w 
and FSSL in eq. (15). (The value of FSSL may be obtained by equating the 
regression rates in eqs. (15) and (20). This point has been discussed in de- 
tail in Section 4, in connection with heterogeneous reactions in composite 

propellant burning. ) Equation (15) predicts the dependence on the initial 

1 

temperature T q as r « 1/ <x [T^/ (T^-T^f . The 5p curve in fig. 3a is 

plotted in fig. 3b for three different values of T q . The recent experimental 

26 

data of Boggs and Zurn are also presented. We note the very close agree- 


ment. 
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IV. TIME- INDEPENDENT BURNING OF COMPOSITE PROPELLANTS 
4. 1 Introduction 

In this section an attempt is made to apply the hypotheses and prin- 
cipal conclusions from the previous sections to the problem of composite 
propellant burning. We continue to consider only those cases in which the 
chemical kinetic rates of combustion reactions in the gas phase are far 
greater than transport and mixing rates; the overall rate-controlling reac- 
tions (viz. fundamental degradation of AP ) occur in the condensed phase. 

The concept of fragment size vaporizing (FSV) is generalized to include 
multicomponent vapor phase/condensed phase equilibrium. That is, the 
contribution of the binder specie s to the vapor pressure is explicitly included. 
In the process, it is required to include the degradation kinetics of the binder 
as well. Thus, the formulation of the problem is "complete. " The influence 
of burning rate catalysts is examined in the light of the same model used for 
uncatalyzed propellants. The increased reaction rate in a surface layer ac- 
counts for both the increased regression rate and the pressure index n. 

The effect of oxidizer particle size is studied in the simplest possible man- 
ner. In the condensed phase, its effect is introduced as an increased surface 
area for the degradation of AP. In the gas phase, the effect is anticipated 
as affecting the spatial extent of the combustion zone, although a quantitative 
formulation has remained elusive at the present time. The regression rate 
is predicted (by a slightly modified form of equation 15) as a function of the 
chamber pressure and the wall temperature. The case of subsurface reac- 
tions only and the case of surface reactions augmenting subsurface reactions 
are treated separately. 
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In order to render the system self-determined, it is required to 
specify the waill temperature. through gas phase energetics and fluid me- 
chanics. For this purpose, the familiar one-dimensional gas phase model 
is invoked. Significant dimensionless groups are evolved and the basic re- 
lation among them is obtained through a formal solution. In the absence of 
a thoroughly satisfactory theory for the gas phase details, a study is made 
of several ways in which the system may be closed. Procedures that have 
been usually successful in determining functional relations in a physical 
problem have limited success here. This is because of the extreme incon- 
sistency in the magnitudes of wall-temperature variation effects in the py- 
rolysis law (strong exponential ) and fluid-dynamic heat transfer (weak loga- 
rithmic). For the purposes of further calculations here, it is seen that the 
assumption of constant wall temperature leads to results sufficiently close 
to reality. 

4. 2 Rate -Controlling Reactions in the Subsurface Region 

We start with the assumption that the rate- controlling reaction (AP 
degradation) occurs in a thin layer on the surface of AP crystals in the solid . 
The propellant material is represented as homogeneous for the purposes of 
heat transfer calculations. That is, a quantity like temperature is assumed 
to have a meaningful and unambiguous interpretation at any plane. The 
process of heat conduction is assumed to be adequately represented by 
Fourier's law with a properly averaged material property of thermal con- 
ductivity. The rest of the analysis will be very similar to that for pure AP 
(Section 3. 2) except for a detail noted below. 

When we consider the differential element dx in the analysis, and 
write the reaction rate in that element, we have to properly consider the 
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actual volume in the differential element in which degradation reactions 
are taking place. As a first approximation, am unimodal distribution of 
spherical AP particles in the solid is assumed. The diameter of AP parti- 
cles is a . If the reactions take place in a layer of thickness t]i on the sur- 
face of AP particles, the volume of AP undergoing degradation reaction 

2 3 

per particle is ira t|i . The physical volume of an AP particle is ira / 6. 

If v is the volume fraction of AP loading in the propellant, the volume 

fraction of AP undergoing degradation is 

6 iIjv 

a ’ 

so that the reaction rate in the layer of thickness dx is modified by this 
factor. We may now absorb this factor into an "effective" pre- exponential 
factor in the reaction term 


B 


- B -(^) • 


( 21 ) 


Recognizing that the rest of the analysis is identical with that for a homo- 
geneous solid, we substitute equation (21) into equation (15) and obtain 


■ {(*£) 


hB* exp(-E/RT m ) ^ j 
(E/RT )X-X 


y 


( 22 ) 


It is worth remembering that the coefficient B is directly propor- 
tional to pressure. The value of the interfacial layer thickness i|f has to be 
specified from considerations external to the analysis. At this stage we 
leave it as a parameter. Given a valid prediction method for the fragment 
size vaporizing (FSV) for the composite propellant, equation (22) may be 
evaluated at any desired experimentally- determined condition. 

Fragment Size Vaporizing. The vapor pressure criterion is directly 
applicable, irrespective of the nature of the vapor species. That is, given 
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the chamber pressure and the wall temperature, the rule in Appendix B may- 
be used to predict the FSV as in other cases. It is only when one asks for the 
value of the wall value of AP fragment size in the burning composite propellant 
that careful interpretation becomes necessary. The statistical mean fragment 
size has to include averaging over the binder species also. 

If m^ gms. of binder (of mean molecular weight FSV^) and rm, gms. 
of AP fragments (of mean molecular weight 117. 5 XFSV 2 ) are present in the 

vapor phase at the wall, the mean molecular weight is 

A m . m_ 

M x - FSV 1 + 117. 5 x FSV 2 ) ‘ (23) 


Under the basic assumption of a plane wall surface (i. e. , uniform consumption 
rate of fuel and oxidizer), the weights m^ and m^ are fixed by the initial mix- 
ture ratio in propellant formulation. 

The assumption of uniform regression rate also furnishes another re- 
lation between FSV^ (fuel) and FSV 2 (AP). The regression rate for the binder 
alone has to equal the regression rate for AP alone (and both have to equal the 

regression rate for the propellant). The regression rate of binder alone has 

1 8 

been the subject of an extensive study in the past ’ and is predicted by eq. (15) if 
the thermochemical constants used in equation (15) are those of the binder ma- 
terial, including the binder fragment size vaporizing at the wall. 

That is. 


r AP r binder (24) 

is the equation that removes the arbitrariness in the values of M and M_ . 

X 6 

For a typical binder (CTPB) and typical numerical values of composite 
propellant characteristics, the predictions are presented in fig. 4. 

4. 3 Joint Rate Control by Subsurface and Surface Reactions 

There are many cases of practical interest where the assumption of 
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subsurface reactions alone is insufficient to account for all of the observa- 
tions. It is suggested that under certain circumstances surface reactions 
in a thin layer augment subsurface reactions. It is found that many of the 
experimental observations can be coherently interpreted. 

The physical model is shown in fig. 1 c. After the material reaches 
the plane SL in its travel from the deep solid region (oo) to the surface, 
the binder and the oxidizer mix thoroughly in a melt layer of thickness L. 
This well-stirred surface layer is at a constant temperature T . In case 
the propellant has burning rate catalysts added in it, the catalyst is taken 
as inactive in the subsurface region (below the plane SL.) and active in the 
surface layer. Such a model is presented on physical grounds. The catalyst 
for its action has to mix well at the microscopic level, a process that is dif- 
ficult to visualize in the solid, but seems natural in the melt layer. 

Equation (22) is directly applicable at the solid-liquid interface 
plane SL. We have to use the value of the AP fragment size FSSL in the 
expression for the non- dimensioned eigenvalue X. However, for a quanti- 
tative determination of the value of FSSL we need another equation relating 
FSSL to the regression rate. We note that the regression rate determined 
by the melt layer reactions (i. e. , without involving any detail below the 
plane SL) is simply given by equation (16). The value of the average num- 
ber of bonds N is given by equation (17). That is, the regression rate is 

L - l 1 - (wv + rek )/ 2 } CB - e*P<-E/* T w >] 

r = 1 r ~ 1 " 1 • (25) 

l FSV " FSSL J 

Obviously, the regression rate determined by the subsurface reactions 
(eq. 22) has to equal the regression rate determined by surface reactions 
(eq. 25). Equating the two regression rates, we solve for the unknowns 
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FSSL and r simultaneously. For the case of catalyzed propellants we need to 
use the appropriate value of the reaction rate expression within the square 
brackets, in eq. (25). 

The numerical solution for a typical case is presented in fig. 5, 
where the thickness of the melt layer is determined. The regression rate 
is plotted as a function of pressure in fig. 6 for the same case. 

Discussion. For a set wall temperature and chamber pressure, 

there are two distinct ways in which the surface degradation contribution 

can be enhanced. An increase in the melt layer thickness L is one, and 

an increase in the reaction rate B* exp(-E/RT ) is another. However, the 

r w 

ultimate effect of either is felt the same way in the regression rate (see 
equation 25). That is, curves of regression rate versus pressure with L 
as a parameter may be considered the same as curves with the surface 
degradation rate as a parameter (possibly through catalysts). 

Figure 6 displays curves where the pressure index n is greater 
than 0. 5 . It is also seen that practically any value of n (within the funda- 
mental range, that is 0 to approximately 0.8) it is possible when one consid- 
ers the possibility of surface reactions. This provides a very direct expla- 
nation of the experimentally observed values of n with catalyzed propellants. 

The variations of the melt layer thickness L with regression rate 
(pressure) will have to be taken into consideration in a more refined analysis. 
We can, however, anticipate the qualitative trends. Since it is known that 
the melt layer thickness decreases with increasing regression rate, the 
regression rate versus pressure curve is likely to be of the shape shown by 
the broken line in fig. 6. We recall that propellants with binders that melt 
readily have been found ^ to exhibit such a decreasing n behavior. 
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4. 4 Gas Phase Details 

All along, it has been assumed that the wall temperature T can 
be prescribed from measurements or from calculations external to con- 
densed phase analysis. In any combustion situation, however, the wall 
temperature is determined by the interaction of the solid phase with the 
general fluid dynamics and energetics of the gas phase. For premixed 
flames (pure AP, and double-base propellants) the full solution is available 
in the literature and will not be pursued here. For non-premixed flames, 
as in the burning of composite propellants, the detailed solution to the gas 
phase depends on the assumptions we make regarding mixing and com- 
bustion. Several considerations are worth our thought before a theoreti- 
cal analysis is undertaken. 

For convenience, the gas phase is analyzed within the framework 

of a one -dimensional analysis. Under the basic assumption of chemical 

kinetic rates being far higher than mixing or diffusion rates, the main 

combustion reactions may be taken as occurring in a region that is standing 
❖ 22 

off a distance X from the propellant surface. Physically, the concept of 
flame standoff distance is as (un)realistic as the concept of a plane wall at 
the surface of the regressing propellant. However, in the analysis, the 
former concept is as useful as the latter. The overall model is depicted 
in fig. 1, which is also found in the current literature. 

Now we turn our attention to the mechanistic details. It is difficult 
to characterize the gas phase above a burning composite as either laminar 
or turbulent. For double-base propellants ( and for composites having 
very small oxidizer particles) there appears to be little room for doubt 
regarding the existence of laminar flow, simply from Reynolds number 
’'See refs. 16, 22, 32. 
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considerations. For normal composite propellants, photographs of the 
combustion region reveal that the gas phase is far from being laminar. 

There are present "spots" which give rise to transport and mixing on a 
scale far larger than molecular. However, the flow field may not merit 
the use of the word "turbulent" in that the origin of the large-scale spottiness 
has little to do with classical Tollmien-Schlichting instability. Moreover, 
and this is an important consideration, it is difficult to conceive of an ex- 
change mechanism which would establish the standard turbulent energy 
spectrum that incorporates proven features such as, for example, Kol- 
mogorov and Heisenberg limits. These thoughts render inapplicable a 
host of useful empirical rules that fluid dynamicists have evolved on flows 
that are truly turbulent. The point to note is that we cannot assume the 
flow to be turbulent simply because it is not laminar. 

In an attempt to avoid these complications, if one examines a single 
crystal of oxidizer and details of gas evolution and flow, the simple laminar 
diffusion flame analysis may be utilized. However, the parameter 

u _ (mean velocity of gas flow)- (flame standoff distance) 

® (coefficient of interdiffusion of fuel and oxidizer gas) 

3 3 

has a value usually less than 5 so that the classical Burke-Schumann 
model without consideration of axial diffusion is of limited validity. For 
the purposes of steady- state heat transfer calculations, it is sufficient to 
work with the non-dimensional flame standoff distance Q = uX / n as a 
parameter in the analysis. A more detailed study , required in the case of 
oscillatory burning, is deferred to Section V. 
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Analysis of Gas Phase Energetics . As already stated, the solution 
depends on the assumptions we make regarding combustion in the gas phase. 
In the absence of a complete understanding of the physical processes, two 
limiting cases are considered. 

Case (i) . The familiar flame-sheet approximation; no combustion until the 
gases have traversed a distance X ' 1 from the surface, and complete com- 
bustion in a short distance after X 1 . 

3 

Case (ii). Uniform combustion at a constant rate m 1 " (gm/cm sec) start- 
ing from the solid-gas interface plane ( ) 

W 

Physically, case (i) is probably a close approximation when the 
"delay" in combustion could be due to one or more of the following reasons. 
An "ignition" temperature has to be exceeded before combustion can com- 
mence. Due to macroscale coarseness of gas phase mixing, a certain dis- 
tance has to elapse before the molecular mixing (needed for chemical reac- 
tions) can be attained. Even if molecular mixing started at the wall, certain 
stoichiometric proportion has to be reached for combustion to initiate and 
sustain itself. 

Case (ii) is probably a satisfactory representation when the mixing 
process controls combustion as in the familiar concept of "well stirred re- 
actor. " Under the basic assumption of molecular diffusive processes con- 
trollingmixing (and hence combustion), the reaction rate may be taken as 
directly proportional to the density of the gases. Since it is known that the 
combustion temperature does not vary appreciably with mean pressure, it is 
easily seen through the state equation that the reaction rate m 1 " is directly 
proportional to the pressure. The molecular mixing rate also depends on the 
interfacial area of contact between the fuel and oxidizer gases and on the ve- 
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locity gradient between them. These inertial processes could be affected by 
the mean velocity of flow. For the present, the interfacial area and velocity 
gradients are taken as independent of velocity (and hence independent of the 
mean pressure through the law r = ap n ) so that the overall reaction rate m'" 
is directly proportional to the first power of mean pressure. 

Case (i) The Flame-Sheet Approximation 

The energy equation applied to gas flow between the wall plane w(-) 


and the burning plane b(+) is 


. d 2 T , dT n 

k + puc -r- = 0 , 

dx QX 


(26a) 


with the boundary conditions 
x = 0 : 

and 

x = x* : 


T = T. 


(26b) 


T = T 


w 


T, - T 
b w 


The solution to the above linear system is written readily as 

T = T • [l - exp(-ux / h )} 

1 - exp(-ux*/n ) ® 

S 

(k is the average thermal diffusivity for the gases, h = (k/pc) ). 

S' s s 

Simple mass continuity gives 

p r = p u . 

H s K g 

This equation may be used to write equation (27) as 

T = T b - f [1 - exp(-rx/H)} , 

with the nondimens ional flame standoff distance 

C = rx^/n 

s 

The temperature gradient at the wall is given by 


(27) 
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/ dT \ _ r_ . v exp(-C) 

\ dx / x ' b~ w' 1 - exp(- C ) 
w 


Defining 

T * (T-T o )/(T w -T o ), 


(28) 


y = rx/ x , 

and 

W = (T, -T )/(T -T ) , 
b w wo 

we write equation (28) as 

-(§) = • < 29 > 

V/ 

If Q (cal/gm) is the heat released by the gases upon complete combustion, 
the overall energy balance between the planes w(-) and b(+) gives 

- k (S) = Qpr . pcr(T b -T w ) . (30) 

w 

Equation (30) may be used to evaluate the flame temperature 
through the equation preceding eq. (9), since the LHS of that equation has to ■ 
equal the LHS of (30). A more convenient procedure is to regard as de- 
termined from fundamental thermodynamics and leave eq. (30) out of con- 
sideration. 

Equating the nondimens ional wall temperature gradient as determined 
by solid phase details (eq. 9) to the same quantity as determined by the gas 
phase details (eq. 29), we write 

C = In {l + 1 + h/FSV } • (31) 

An examination of the fundamental processes from several viewpoints con- 
sistently leads to the same groups, W , (1 + h/FSV), and W/(l + h/FSV) . It 
is suggested that these four groups be regarded as fundamental similarity 
parameters in propellant burning, in the same class as the groups 9 and 

Si 
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X in solid phase analyses. 

To have a feel for the variations in £ with regression rate and pres- 
sure, fig. 7 was prepared. The curves are computed through eqs. (15) and 
(31) for the case of a composite propellant with typical, assumed values of 
the various parameters. The binder contribution to the vapor pressure was 
neglected since such details are thought unimportant in determining the £ 
trends. 

It is interesting to note that (at constant T^ or an increasing T^) 
the physical flame standoff distance x* is predicted to decrease with in- 
creasing pressure. (This is to be recognized from fig. 7, after noting that 
u / k depends on pressure like p n where n is the pressure index in the law 
r a p n . ) This decrease in x* with pressure increase is directly supported 

5 

by experiments . While it is true that such a decrease in x* is also pre- 
dicted by most other contemporary gas -phase theories, the result depends on 
the gas-phase chemical kinetic reaction rate in those theories. In terms of a 
simple picture, the chemical reactions can be completed in a shorter distance 
at higher reaction rates (at higher pressures ) , whereas, in the present case 
we see that the result follows even when the gas-phase chemical reaction 
details are completely ignored from the analysis. 

It may be seen from figs. 3-7 that the assumption of constant T 
leads to predictions that are sufficiently close to experimental data. Further 
discussions on the subject of wall temperature determination through gas- 
phase details may be found in Section V. 4 . 

Case (ii) Uniform Combustion 

The solution here requires, for numerical computations, the value of 
the reaction rate m'" , which has been taken from the work of Culick and 
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34 

Dehority . Also, in order to preserve the similarity to that work for easy 
comparison, the coordinate system is now taken as, distance x measured 
positive from the wall plane in the streamwise sense (see fig. Id). 

The energy equation and boundary conditions are: 

. d 2 T dT A . tlI 

k T - 0 c u — — = - Qm 1 " 

g ^ g P dx 


dx 

x = 0 : T = T 


w 


(■^1) = p uc(T -T ) + -=JL p 

\ dx / *g w o FS V K 

ttr 0 


u 


w 


X = X* : T = T, 


g 


<&>. 


J 


(32) 


The definitions 


T - (T-T 0 )/(T w -T 0 ) 


y = (ux)/h 


g 


m = p u 
g 


V * 111 

A 2 = g Qm 

A g It (T -T ) 
0 me w o 

P 


C = (ux*)/k 


g 


transform the system (32) to 

d 2 T dT _ .2 


g 


y = 0 : t 


1 


il = 1 +-JL. 

dy FS V 

y=C: T=l + W 
dT 


-\ 


(33) 




dy 


= 0 (inapplicable if A = 0) 


( 34 ) 
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The symbol C is not the same as used by Culick and Dehority ; C i n ref. 34 

would be exp(y) in the present notation. 

The solution to the linear system (34) is written readily as 

t = 1 + A 2 {y + exp(-0 exp(y-C)} . (35) 

S 

Recalling the earlier definition, 

W s (T u - T )/ (T -T ) , 
b w w o ’ 

eq. (35) may be made to yield 

O = < 1 + - exp(-C)} (36) 

7 w 6 

and 

T x=x * = U+W) = 1 + Ag(C + exp(-C) - 1} . (37) 

In principle, the problem is completely solved. For any prescribed 

pressure (and hence m'" ) a trial value of wall temperature T w is assumed. 

The mass burning rate, m , is computed through eq. (22) (or eq. 2 5, if 

2 

that case is being considered). The value of A is computed through eq. (33). 

§ 

The nondimens ional distance Q is computed from eq. (36). The heat release 
parameter W is now computed from eq. (37). However, the parameter W 
can also be computed from the overall energy balance equation (30) which is 
written in nondimens ional form (after the use of earlier definitions) as 

(1 = — (y 9 :t t - w • < 38 > 

p w o 

The iterations are repeated until satisfactory agreement is found between 
the two independent values of W. The results are presented in fig. (4b) for 
a case where the binder details have been ignored in the condensed phase 
treatment. It is a straightforward procedure to include them (as done in 
fig. 4a), but the details were thought to be unimportant in the present context 
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where the aim is the elucidation of gas phase details; besides, the values of 
Q (or ) are not readily available even for the simplest of specific pro- 
pellant compositions. 

In fig. 4b the burning rate pressure index is seen to be 0. 46 . In 
comparison with predictions based on the assumption of constant wall tem- 
perature (see fig. 4a), the general behavior is seen to be not appreciably 
different quantitatively, or even qualitatively. 

The variations in wall temperature and the flame standoff distance 
(strictly speaking, the plane of onset of zero temperature gradient) are 
plotted as a function of pressure in fig. 7b. Again, in comparison with pre- 
dictions based on the assumption of constant wall temperature (see fig. 7a), 
the general behavior is seen not to be appreciably different. 

It would appear that the actual model for the gas phase processes is 
not crucially important in determining the general trends so long as the 
overall rate controlling reactions occur in the condensed phase. Such a 
conclusion was anticipated earlier (Section II. 3) purely on physical grounds. 
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V. OSCILLATORY BURNING OF COMPOSITE PROPELLANTS 
5. 1 Introduction 

The role of the propellant in pres sure -coupled instability in rocket 
motors is thought, at the present time, to be well represented by the response 
function, defined as 



with the familiar nomenclature. It has been shown that the numerous the- 
oretical derivations of the response function are essentially identical, since 
the fundamental assumptions in the various analyses are identical. 

On the experimental side, some doubts exist regarding proper pro- 
cessing of the measured variables to yield the real part of the response 
function. Direct measurements of the fluctuating regression rate are dif- 
ficult, and the several different techniques that are in use have not proved 
conclusive. Secondly, measurements at low frequencies involve large com- 
bustion chambers so that the question of thermal inertia and equilibration 
has introduced some doubts regarding the acceptability of data in the low fre- 
quency regime. Nevertheless, some of the general trends in the experimental 
response functions are thought to be well understood. The response function 
is usually found to be dependent on the mean chamber pressure, contrary to 
theoretical predictions. In some experiments, the response function seems 
to assume very large values at low frequencies, while theoretical predictions, 
and even physical arguments, suggest that the response function must tend to 
the limit of the steady- state pressure index n, in the limit as frequency tends 
to zero. Reference 37 presents an excellent introductory treatment of the 
general problem of combustion instability. 
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Based on the present model of composite propellant combustion, 
theoretical derivations of the response function are undertaken in this section. 
The cases of subsurface reactions only and surface reactions augmenting 
subsurface reactions are treated separately, since a cursory examination of 
the voluminous experimental data on response functions indicates that such a 
classification of propellant burning may be an important one. The analytical 
solutions are derived in a manner very similar to the steady-state analyses 
in the previous sections. The theoretical response function is found to be a 
function of the mean chamber pressure. Within the framework of the adia- 
batic approximation for flame temperature fluctuations, the theoretical re- 
sponse functions (for the case of subsurface reactions only) are seen to ap- 
proach large values at low frequencies. However, it is recognized that 
ultimately the adiabatic variations have to give way to isothermal ones, thus 
removing the inconsistency in the low frequency limit (at zero frequency). 

In a second attempt to postulate a physically valid model for the gas 
phase processes under oscillatory combustion, uniform combustion is con- 
sidered in the gas phase. The general trends in the response function are not 
found to be very different from the case of adiabatic approximation, although 
the mean pressure dependence is now found to be much weaker. 

The theoretical response functions derived in the present study are 
qualitatively discussed at the end of Section 5. 3. It would aid the reader to 
study the figures 8, 8b, 9, 10, and 10b after going over that discussion. 



-41- 


5. 2 Subsurface Reactions with No Surface Reaction . 

The real problem in the solution of the time dependent case, it turns 
out, is the specification of the proper boundary conditions on temperature 
gradients. The rest of the procedure is very similar to the time independent 
case discussed in Section III and needs no elaboration here. 

We recall that the temperature change across the reaction zone is of 
the order of reciprocal activation energy parameter, i. e. , ~ 1/9 . In terms 

3 . 

of physical distance, x, this reaction zone has a thickness ~ h/ r0 . Thus, 
we see that, so long as we restrict our attention to small amplitude (~ 1/fi ) 

cl 

fluctuations, the effect of temperature variation across the reaction zone 
is a second order quantity and need not be included in a first order analysis. 
That is, if we visualize a (hypothetical) freezing layer below which chemical 
reaction rates are negligible compared to other processes, we can neglect 
the effects of mean temperature variations between the wall and the freezing 
layer on fluctuations of temperature. 

It will be shown (see later) that phase differences in variables may be 
neglected across the reaction zone except at very high frequencies. This 
fact can be exploited to write the exact expression for the wall boundary con- 
dition. 

We shall first use a coordinate frame of reference held fixed at the 
mean position of the fluctuating wall (or freezing layer ). We obtain the outer 
solution, neglecting the reaction term. We then use this exact solution to 
write the temperature-gradient boundary conditions both at the freezing 
layer (the interface between the "inner" and "outer" regions) and on the wall. 
Having obtained the boundary conditions, we solve for the time -dependent 
regression rate in a manner analogous to the time -independent case. That 
is, we solve the full equation, including explicitly the Arrhenius reaction 




-42- 


rate term in the inner region. 

Outer Region. In the outer region, where the reaction rate term is 
exponentially small, we may write equation (1) as 



, - 3 T 

+ P cr "EE 




(40) 


Taking r = r + r' , and defining 

T = (T- T )/(T -T ) = (T-T )/(T -T ) + T/(T -T ) = 7 + T 1 , 
owo o w o 'wo 

y = xr/x. , 



r 


and assuming harmonic fluctuations in t', we separate equation (40) into 


and 


*4 + 21 = o 

17 dy 


24 + *ll . in + = o 

17 d Y 


The solutions to equations (41) and (42) are, respectively, 

T = exp(-y) 

and 


T' = T^,exp(-X 1 y) 


where 


\ l = i(l + a. l + ibj ) 

a x = [(1 + l6C?)Kl¥/ JT 

b 1 = [(1 + 16Q 2 )^- ll^/ V2 


(41) 


(42) 


(43) 

(44) 



-43- 


Since the position of the wall plane (and also of the freezing plane) 
is fluctuating in time, the coordinate there is x^ = JVdt , which may be 
written v 

1 T>l \ 

■ (45) 


w 


. 1 r* 

-1 n - 

r 


Expanding the temperature in Taylor's series around the freezing 

2 

layer, and using equations (43), (44), and (45), it is a matter of algebra 
to arrive at 


-( k iD = r p[Vf + 


T,-T 

for' 


1 r 

which may be written in terms of nondimensional variables as 

1 r' 


-( 


8t y _ \ ~ * 

3y/ f l T w X. — 
1 f 1 r 


(46) 


Equation (46) specifies the outer boundary condition at the inner edge. The 
inner solution, which will be derived shortly, will have to match (in its outer 
limit) with equation (46). 

In order to examine the order of various terms, we write equation (1) 


in non-dimensional form as 


9t 


9t,8t 

9y 8 (art) 


= AhN exp(-E/ft T) . 


(47) 


In the linearization of the fluctuating part of the reaction term, there 

is the inherent assumption that 0 XT' ~ e , which implies that t' ~ . 

aw w 

Now, examining the terms in equation (47) in the inner region, we see that 

9t 


n 


9(u>t) 


is of the order Qe . Hence, the term may be neglected from the equation 
except when 0~ 1/e or higher. Excluding such high frequencies from con- 
sideration, we write equa tion ( 1 ) as 

The symbol e denotes any (arbitrary) first-order quantity. Examples, 1/0 X, 
P'/P, r'/r , etc. denotes, thus, any second-order quantity. 
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k — — 2 + prc 4^ = DpNB exp(-E/RT) . 
dx 


(48) 


Physically, this tells us that the reaction zone is of such a small 
thickness that effectively the phases of various quantities of interest are 
left unaffected across this layer. Viewed in another light, the length scale 
associated with the oscillations (~r/u)) is far larger than the length scale 
associated with the degradation reaction ( — k/ 0 Xr ) . These observations 

cl 

introduce remarkable simplifications in the analysis to follow. 

First, we may use the convenient wall-fixed coordinate system and 
closely follow the (already derived) steady-state solution, recognizing that 
the variables have an implicit time dependence and that new boundary con- 
ditions are to be used. Second, we write the wall boundary condition in 
analogy with equation (46) [see figure le] as 




L — fi + h N 
, - \ F Wv) 


\dy/ ^ 1 T w + XT 

w 1 r 


Lastly, we may substitute in the inner region, 

d 1 d 


(49) 


dy ( 1 + r ' / r ) ^ 

where 

y = rx/ K . 

(We recall that r = r + r' . ) 

Boundary Conditions on Temperature Gradients. 
Freezing layer ( )^ : 

(#P f = (f ) f + O' ^-(VwV^) 


(50) 


or 


(T) = 

’ j * e 


dy 


fl+r'/r } f lr r 
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The variable £ ~ + is written 

dy 

£° = = % r + £' = 1 + T 

b b£ b f b f w 


' - 1 _ X T' - — (l - y!-) = T' (1-X ) + — (l- -rM. (51) 

W 1 W — \ A. , / W 1 — \ A-,/ 

r 1 r 1 


Wall ( ) : 

w 


<£)_ • Q_ + (&)_ = 


rdT 


3t 1 


1- 


dy/ V 3 y / ~ FSV 

w w w 


-Vw " ^ 0 + Wv)xr 

r 1 


or 


(7) -77V =[-( lt ®)-v;T = ('^)](‘-r) 

dywl+r' / r 7 w , 1 r ' r 

1 

i 

= -( 1+ wv)- x rw-i7^( 1 + i5v) + ^ ( 1+ ilv) • <52) 

' 1 r r 


Now we have 

| =1 + %' = - + T' (1-X. ) + — (l + -s—Vl - 77 -) . 

*w ^w w FSV w 1 — \ FSV/\ A., / 

r , 1 

j 

It is convenient to employ the shorthand notation 


w 


— — + A + — — 
FSV 5 FSV 




1 r 


(53) 


(54) 


which is nothing more than a substitution; of equation (51) into equation (53). 

Now we have all the information necessary to solve the time -dependent 
regression rate r . 

At the very start we observe the satisfactory feature that the boundary 
conditions, eqs. (51) and (54), tend to the proper steady-state limit as the 
frequency -* 0 (i. e. , as X^ -*1). In all^of the analysis to follow, it is worth 
remembering that the variable T and the parameter D/ c have been normal- 
ized with respect to the steady state quantity, (T -T q ) . 

The Solution. In the inner region, where the gradients get very 
large, we write equation (12) as 


• V* 
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pp'+p = (h+p+T ) A exp[-0 { 1+x(1-t )}J . 


(55) 


Using the small parameter e = 1/ A X to stretch the variables in the inner 

cL 

region, and defining the inner variables 


■n = 


_ 1 - T 


and 


| = p + T 


X = e A exp(- 9 ) , 


we write equation (55) as 

-(| - 1 + e r i)£' = (h+^)^exp(-r)) 
where a prime denotes d/dr| . 

Making the expansions 


and 


£ = £ 0 + ®£ 1 + e l 2 + . . . 


r = r + e , x 1 + e %+ ... 

o 1 2 


(56) 


and substituting them into eq. (56), we get, to the lowest order, after dropping the 
subscript o, e' - £'£ = (h+£ )^exp(-ri) . (57) 

To clearly display the integration we write equation (57) as 

* X exp(-Ti)dr] , 

which integrates to 

An(h+£ )- [£-h 4,n(£+h)] = -X exp(-ri) + [constant! , 


i. e. , 


(h+1 )£n(h+£ )-£ = -X e -r ' + [constant] 


(58) 


( 59 ) 


Kaplun's matching theorem requires 

t t° 

| -* | as p -* oo , 


yielding 
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(h+1) -tn(h+£ 0 )-£° = [constant] , 

f. o 

(h+l)<tn{h(l + -£° = [constant! , 

i. e. , 

(h+l)tnh + (h+l)-tn(l + £°/h) - = [constant] . 

Recognizing that £° ~ e , we approximate the logarithmic terms and 
write 

(h+l)-tnh + (h+l)£°/h - £° = [constant] 

(h+l)-fnh + £°/h = [constant] . (60) 

Equation (60) may be substituted into equation (59) to yield 
(h+lRn(h+£')-| = -r e -T1 +(h+lRnh+ £°/h 
or 

X e = (h+1) ■tn + £ + £ ° / h . 


On the wall. 


T| = 0 


and 


e = e 


w 


hence. 


X = (h+1) <ln 


h+e 


+ i + e°/h 


w 


(61) 


r 

FSV ' 

h r ' (l- 1 ^ 

r h+i 

h+1 

h "I 

L 

FSV- 1 . 

FSV r V 

L h(FSV-l) 

FSV-1 

FSV -1 


After considerable algebra, detailed fully in Appendix A, we write 
equation (61) as 

5V.-1 _h__£l(l 

7-1 J FSV r V 

" T w (X l -1) h(FSV-l) * (62) 

Fluctuations inwall temperature, and chamber pressure fluctuations, 
produce fluctuations in the value of FSY. . Its effects are second order in the 
third and the fourth terms in equation (62) and have been neglected. It is 


h+1 



shown in Appendix B that FSV can be satisfactorily represented as 

FSV = a P _(3 • exp(6T w ) 

where a, (3, and 6 are constants. 

After some algebra detailed in Appendix A, we write equation (62) 


as 


X = (h+l)-tn 


FSV 


+C, 


El 


FSV-1 FSV 


r-( 1 -ir) c 3- T w' c 4 +c i (X r 1)5 
r 1 


(63) 


where 


C 1 “ 


C 2 5 


h+1 


h(FSV-l) 
h+1 h 


FSV-1 FSV 
C 3 * C 1 + C 2 


C 4 - C 2 6T w X 


C 5 " C 2 P 


We recall (equation 14) that under steady conditions 
X = (h+lRn[^V]- h 


FSV ’ 

which can be used to normalize the unsteady \ (equation 63). 

X 

- 1 - nr T „, T T- r- - — V 1 - 1 , 

p r 1 


X . 4 , C 5 2L r ' A M 3 t\ n 1 

X = 1 - -JT T w + t - - ~ V 1 - rbr t • 


( 64 ) 


It is worthwhile to recall the original definitions, 
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X = 


x B(1 + p'/p)‘ exp(- 0 ) * (1 + 0 XT' ) 
a a w 

9 U+Xt' )fx( 1+t' (l-X)}](r’+r l ) 2 
aw w 


The expressions for X and X give, through equation (64), 
ft 


i - c 5 /x 




P /p 2 - (l-l/XjC./X - — 2L. JT X+C./X-1+(X -1)C./X] 

1 3 r'/r a 4 11 


(65) 


Equation (65) is the desired expression for the response function. Given a 
relation between and r'/r , the problem may be considered "solved. 11 

A Relation Between t' and r'/r . A valid procedure is needed for 
specifying the relation between and r'/r , external to the solid phase 
analysis. For a propellant regressing under the influence of its own gas 
phase reactions, it is logical to invoke the gas phase details to supply the 
relation between and r'/r. However, we recall that a thorough and 
complete solution to the gas phase has not been obtained yet even for the 
time -independent case. For the present, we shall be content with choosing 
any one of the several criteria that were evolved in an attempt to eliminate 
the arbitrariness in time-independent burning. The data in fig. 4a suggest 
that constant (mean) T w is an adequate assumption, reasonably representa- 
tive of many experimented measurements. From dimensional considerations 
we expect that C, = rx*/n does not vary too much from its steady state value 
during unsteady combustion. It is assumed in this work that C is identically 
equal to its value during time -independent burning at the same mean pressure. 
The variations in the value of C with mean chamber pressure have, how- 
ever, not been neglected. 

The temperature and pressure fluctuations have to match the "far 
field" values just beyond the "flame sheet. " Gas phase acoustic processes 
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in the main chamber are usually fast enough to be isentropic (but not so fast 
as to introduce effects of thermodynamic nonequilibrium! ). A relation be- 
tween the temperature and pressure fluctuations follows immediately as 

^ = 3LdEl . (66 ) 

A " V 

Some caution is necessary in the use of this seemingly universal equa- 
tion (66). When we consider the low frequency limit, (JU -* 0 , the above equa- 
tion is obviously not representative of reality. It is known that flame tem- 
peratures vary little with mean pressure in time -independent burning. For 
a very slow change, the process is more likely to be isothermal than adia- 
batic. A convenient way of handling this effect mathematically is to let y -*1, 
in eq. (66), as uu-*0. In an ideal analysis, the index z in the expression 

T' /T = — — p'/p, will be written as a function of frequency such that z -*1 
z 

as uu -* 0 and z-'y, sufficiently far away from uu=0. It is found that a very 

3 

similar argument has been put forward by other researchers , who have ac- 
tually presented T'/T = fn(p'/p, uu). However, both the gas-phase physical 

model they use and the analysis are thought to be too involved to be consistent 

3 

with the present work, if we think of using the available expression here. A 
simpler procedure for handling such effect is being considered at the present 
time. 

The debatable assumption of isentropic fluctuations in the gas phase 
may be dispensed with if some other criterion can be invoked to evaluate the 
ratio T^/(r'/r). In the calculations of time-independent regression rates 
(see IV. 4), the model of uniform combustion in the gas phase was invoked to 
dispense with assumptions regarding the value of wall temperature. As a 
generalization from the time-independent case, the same model of uniform 
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combustion in the gas phase can be used in the time -dependent case also. 

We thus have two different models for the gas phase (adiabatic and uniform 
combustion), both of which can be used (independently) with the same con- 
densed phase model. 

(a) Adiabatic Fluctuations in the Gas Phase 

Except for the small intricacy at very low frequencies, eq. (66) is to 

3 6 

be used in the analysis. As discussed by Krier, et al. , many experiments 
are available where the isentropic assumption is valid. Hence, the present 
analysis finds use even without any refinement. 

Under the basic assumption of quasi-steady processes in the gas 
phase, we recall that the value of the wall temperature gradient is 


-(~) 

dy w 


T - T 
b w 

T -T 
w o 


exp(C )- 1 



T -T T -T 
wo wo 


exp(C) - 1 


T 1 

w 

T - T 
w o 


■pi T T T* 

W + — - • — . — - — — at_ 

T, T-T T T-T 

b w o w w o 

exp(C) - 1 


(67) 


Recalling all of the assumptions and nomenclature, it is a simple exercise 
in algebra to write equation (67) as 


/±L\ _ (i + 1 + Izl 2 W ±1 £l . Ik ~| 

_ \ l + fsv'L y xw - w J 

dy w p 


( 68 ) 


Equating equation (68) to the same temperature gradient evaluated through 

solid phase analysis (eq. 52), we have 

t' 

. v + 


-1 X W+ 1 p' w 
— X W - ' W 


1 + 


FSV 




(69) 
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It is a matter of convenience to work with 

h 


C 7 (1 + FSV 


)/2 


^ = y-1 XW + 1 

C 8 y X 

R20 ■ (l ' ) 


which simplify equation (69) to read 


t' 

w 


r’/r 


O 1 

w i - R20 



where R = (r'/r )/ (p'/*p ) . 


(70) 


The value of T^/(r'/r) predicted by equation (70) may be substituted into 
equation (65) to obtain a formal solution, i. e. , an expression for the 
response function, R. For convenience again, we employ short hand no- 
tations : 


R6 = 
R41 = 

R42 = 


(e a x- 1 + c 4 /\ + (X 1 -i)c 1 /\ } 

C 8 /W 

JL +ii_ 

w 2C ? 

R20 

W 2C ? 


R43 

R44 


1 - 2X R2 ° + 

R6 • R41 
z 


R42 • R6 
2 


It is to be recognized that the group (l-C^/X)/2 is the pressure index n 

for time-independent burning (under the assumption of constant mean T ) . 

w 
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Finally, the normalized response function is written 



n 


1 + R44/n 
R43 


( 71 ) 


Discussion 

Before an attempt is made to compare the theoretical curves in 
fig. 8 with experimental data, the validity of the theory must be clearly 
understood. 

The curves are not valid at low frequencies because the index y 
(in T' / T = ■ - p' /p ) approaches unity and does not retain the value of 1.25 

that was used in preparing the graphs. The effect of this shift to isothermal 
changes at low frequencies is to converge all the curves to the ordinate 
unity as Q — 0 . 

The curves are not valid at high frequencies ( Q k, 1 0 ). In the singu- 
lar perturbation analysis we excluded such high frequencies from considera- 
tion. 

There is a third limitation which may not be readily obvious. The 
curves are not valid at high pressures, pr to be more precise, when the 
length scale k/ r becomes smaller than the mean size of the oxidizer parti- 
cles in the propellant. The length n/r is the effective thermal depth in the 
solid. When this depth gets comparable to, or smaller than, the mean size 
of the oxidizer particles, the material cannot be considered homogeneous 
any more. When the assumption of homogeneous solid breaks down, so 
does our analysis. 

7 10 

There are several reported experiments ’ that closely follow the 
trends predicted in fig. 8. The numerical values of the real part of the re- 
sponse function measured in experiments are close to those predicted in 
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fig. 8. The frequency scale appears to be shifted slightly toward the higher 

-3 2 

side in the experiments if we use a value of 1. 1 X 10 cm /sec for the 
thermal diffusivity h, as we have been doing all along. A different value 
(but not unrealistically different) for the thermal diffusivity would remove 
this frequency position discrepancy. 

We regard with satisfaction the existence of experimental data that 
agree with the theoretical predictions. We recognize that the ultimate aim 
is to predict the outcome of any experiment. 

(b) Uniform Combustion in the Gas Phase 

Under the basic assumption of quasi- steady processes in the gas phase, 
the solution to the gas phase energy equation is valid during time -dependent 
burning also. 

The wall temperature gradient was seen to be (Sec. IV. 4) 


• a n = . 


w 


Vdy / 


w 


* g 2 »-e-C) = 


p 2 r^C 2 (T -T ) 

r p wo 


(l-e" C ) 


From ref. 31 it may be seen that the mass consumption rate of gaseous re- 
actants can be represented as 1 

m'" = A • P (72) 

where P is the mean pressure and A is an appropriate constant. 

Defining the nondimens ional mean heat release rate as 


* (A) 


CL 


S c(T -T ) 
w o 


A -P n -C. 

-2 (1_e > 

P r 
r s 

i 


(73) 


the wall temperature gradient during time -dependent burning is written 


p = -S(l + P'/P)(l - 2r'/r) . 

w 


( 74 ) 
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In the above equation, the assumption (as first introduced in subsection (a) 

before) is made that the nondimensional reaction distance Q can be taken as 

constant during unsteady combustion. The variations in C with mean pres- 

\ 

sure have not been neglected. 

It is easily verified that the nondimensional S defined in eq. (73) 
has to equal the steady- state nondimensional wall temperature gradient 
(1 + h/FSV). Although it is a very straightforward procedure to consider 
either, the value of FSV used here is the full time -dependent value and not 
just the steady-state value. This point needs amplification. The heat re- 
lease in the gas phase combustion is not strictly a constant; it depends 
on the initial endothermicity of the propellant vapors. A large value of FSV 
indicates a large number of bonds left unbroken in the vapor molecules and 
the (necessary) scission of the bonds acts like a heat sink on some of the 
energy release in the gas phase due to combustion. Hence, the net , , 
which is what is implied in the formulation here, actually depends on FSV. 
Thus, during unsteady burning, the fluctuations in the value of FSV produce 
fluctuations in the value of . We may- absorb such fluctuations in into 
the (1 + h/FSV) term and leave as a constant in the definition of S (eq. 73) 
Purely from an analysis of the condensed phase, we had the wall tem- 
perature gradient (eq. 52 ) as 

P» = -( 1+ ^V)- x L T w-( ! + Flv)3 0-r-) • 

r 1 

Equating the two separate evaluations of the wall temperature gradient, we 


get 



See Section VI (iv) for a detailed discussion. 
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or 



( 76 ) 


Division of eq. (76) throughout by r'/r , and substitution of the identity 

ft e (r'/7)/(P'/P) yields 


or 



(77) 


so that the problem is solved, in principle. Using the result (77) in eq. (65) 
and employing shorthand notations for convenience, it is a matter of algebra 
to arrive at 


U1 + U2 
~ U4 + U5 


(78) 


with 


U1 = 1 - c 5 /x 



U3 - 1 + l/X. > 

U4 = 2 - R20 • C 3 /X 
U5 = U2 • U3 


(79) 


I 
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5. 3 Surface Reactions Augmenting Subsurface Reactions 

In close analogy with the steady burning case (see 4. 3), the case 
of time -dependent combustion is treated here. The physical picture is the 
same as in fig. lb. The only difference from the steady burning case is 
that the fluctuating part of the regression rates in the subsurface region and 
in the surface melt layer need not be equal to each other, as the steady part 
of the regression rate has to be. The difference between the fluctuating 
parts of the regression rate in the surface and subsurface regions shows 
up as fluctuations in the thickness of the melt layer. 

To anticipate the analysis, the developments in Section 5. 2 are 
carried over as exact representation of the propellant below the plane SL. 
The fragment size is thus FSSL and not FSV ; also, the fluctuations in 
FSSL are not given by a perturbation of the vapor pressure rule B. 2, but 
have to be determined through proper matching with the melt layer. The 
surface contribution in the melt layer is a linear perturbation of the ex- 
pression for the regression rate (eq. 25). The resulting equations are 
solved simultaneously to obtain the response function. The relation be- 
tween the wall temperature fluctuation and the pressure fluctuation has to 
be supplied from gas phase details, as in Section 5. 2. 

In the following analysis, the steady and fluctuating parts are 


written 

surface melt layer thickness 

fragment size vaporizing 

fragment size at the solid- 
liquid interface 

fluctuating part of the regression 
rate due to subsurface re- 
actions only 


L = L + L' I 


FSV = FSV + FSV' 


FSSL = FSSL + FSSL 1 f 


( 80 ) 


r' 


ss 


J 
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It is an exercise in algebra to write equation (62) (" applied below 


the plane SL ] as 


ss 


{2 - 


CL3 


(1 ~ T~ )} = 


CL2 FSSL 1 


FSSL 


T 1 

w 


{(X r i) 


CL1 

~x~ 


■p I 

B X- 1 } + — 
a P 


(81) 


In the above expression the shorthand notations are: 


CL1 

h + 1 


h(FSSL - 1) 


CL2 

h + 1 

h 

“ 



FSSL - 1 

FSSL 

CL3 

s CL1 + CL2 



x „ (h+1 ) in { FSSL — ! h_ 

FSSL - 1 FSSL 

It should be remembered that X^ is the complex root of the unsteady equa- 
tion (see 5.2) and that the above shorthand notations are all steady-state 
quantities only. 

It is again an exercise in algebra to write the unsteady part of 
equation (25) as 


r^ 

r 


FSV' z3 FSSL' 
FSV FSSL 




+ 8 xt 

a 


! 

w 


with 


Z2 = FSV ■ FSS L ( 1 - 1 / FSSL) FSV 

2 • FfTv • F^SSL - FSV - FSSL FSSL - F£fV 

and 

Z3 3 FSV- FSSL (1 - 1/FSV) _ FSSL 
2 • FSV- FSSL- FSV- FSSL FSSL- FSV 


(82) 


both of which are shorthand notations for steady-state quantities. 
Mass conservation in the melt layer gives 
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, _ , _ dL' 

r ss r dt 


( 83 ) 


Assuming harmonic variations in L.' , equation (83) is written in non- 
dimensional form as 


ss 


. Q L'r L w , n L’ rL 
= i — = *-,(*-, - U 


4 Lh 


T 1 


L H 


or 


L/ 

L 


r r' , 


1 


- - J ' 

rL r r 11 

The fluctuating part of FSV may be written through 

FSV = n.P - * 3 exp(6T ) , 


( 84 ) 


w 


as 


FSV' 

FSV 


= C.T' _ C c P'/P 
4 w 5 


(85) 


Substituting the values of FSV'/FSV and FSSL'/FSSL from equa- 
tions (85 ) and (81 ) into equation (82), we get 


^ = Z 2 <C 4 V C 5 P,/P)+ CE2' 

r 


( 86 ) 


T i pi 

+ — + — + e xt' 

L -P a W 


During the discussion on steady-state combustion, the quantity L was sup- 
plied from external considerations. Strictly speaking, it should be determined 
(from a knowledge of the melting behavior of the propellant) as a function of 
FSSL and FSV. in the unsteady case, the fluctuating part L' thus remains 
unknown. We can, however, consider £wo limiting cases to obtain a closed 
form solution. 

(i) The melt layer thickness is constant even in unsteady combustion 
as in steady combustion. Now, = r' and L 1 = 0 . 
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(ii) The subsurface regression rate fluctuation vanishes, i. e. , 

r' = 0 . Now, all of the fluctuation in the regression rate is 
s s 

entirely due to the melt layer. 


The physically real solution exists between these two limiting cases, 

probably closer to the second criterion. 

Both the solutions will be discussed here. 

(i) Case of L 1 = 0 (r 1 = r'); (a) Adiabatic Fluctuations. 

s- s s 

It is easily seen that the wall temperature fluctuation is related to 
the other fluctuations through equation (70). That is, the flame temperature 
fluctuations affect the entire energy balance exactly as discussed in Section 
5. 2 Equation (86) is written 


w 


where 


r'/r 

\ 


Z7 Z5 1 
Z4 Z4 ft 


(87) 


24 S Z2 - C 4- c!T2 Z3 R X l- l,£ T i+ V- 1 ]+ 0 a* 


Z5 s Z2- C r + 


5 CL2 


Z3 - 1 


zts ch z3 ^ 2 -• £ r S(1 -x7» 


Z7 = 1 - Z6 


( 88 ) 


Equating the values of t 1 / (r'/r) in equations (70) and (87), we have 

w 


vri - RZ0 

~T~r- 


Z7 Z5 1 
Z4 Z4 ft ’ 


W 2C, 


yielding, through the definitions. 



-61 


R30 

R31 

R32 
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_ W 


R20 


ZL(± + _^_N 

Z4 \W + 2C„ / 


Z5 W 
Z4 C 
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7 ' C 8 
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(J- + _J_) 
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1 - R32 
R30 + R31 






( 89 ) 


(90) 


(i) Case of L 1 = 0 (r|._ = r 1 ); (b) Uniform Combustion. 

S S 1 ’ ' ■ 1 

From the earlier discussion on the model of uniform combustion in 

the gas phase (Section 5.2(b)) we have an expression for t' /(r'/r)in equa- 

w 

tion (77). Equating the two independent expressions obtained for the same 
quantity (equations 77 and 87), we get 

H = <* + f#v>[]r U3lA i • (,1) 


(92) 


Hence, the response function is 


R = 


(1 + 


(1 + 


FSV 


Z5 

Z4 


FSV 


+ Z7 
x i 


Z4 


Case(ii) r '^ /r = 0 and L'/L 4 0 ; (a) Adiabatic Fluctuations. 

- -1 . S S *' ' ‘ ' ~~ . ■ " ■ I I I. .. I 


Substituting the values of FSV'/FSV and L'/L into equation (82 ), 


we get 


7 = Z2[C 4 T w- C 5f- ] - Z 3 { T w r( V 1)£ r i;0 a’<- 1 ] + ? : }|T 
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(93) 
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Defining 


CT 1 \ 

R531 = Z2* C 4 -Z3{(X r l)-^i + e a X-l}~ +6 & x 


R532 s Z2* C 5 -Z 3 + 1 

(■=) 

R ” 3 E 1 " V \ r , - TT 

we write equation (93) as 


R533 • — = R531 t 1 - R532 — 

— W r= 

r P 


( 94 ) 


(95) 


From considerations of gas phase details only we can write (see Section 
5. 2 ) the wall temperature gradient as 


- t 

' 'll T 


1 + c 


dy 'W 


8 p W 


w }(i + — ^ 


FSV 
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(96) 


Recalling that r 1 is zero, we may write the temperature gradient at the 
s s 


plane SL as a simplified version of equation (52): 
/ dT ■’N /, . h 


- (^) - o 3 

dy SL 


) + X i T ’ 

FSSL ' I 1 


w 


(97) 


The difference between the heat transfer rates across the two planes SL 
and w has to be the sink term due to chemical reactions in the isothermal 
layer on the surface. In terms of nondimensional variables, the sink term 
is 


FSV 


_1 \ 

9 


FSSL 

i 

which is valid regardless of whether the processes are steady or unsteady. 
This is both because we are considering a region of no phase lag and the 
normalization of the variable x into y has the full time- dependent regression 


rate included in it. 
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Equating the difference between the equations (96) and (97) to the 
above expression, we get 


t' 

w 


{ 


FfTv 

~w~ 





1 + 


FSV 


>• 


(98) 


Using the convenient notation 

C„(l + h/FSV) 
R534 s 8 


we combine equations (95) and (98) to write the response function as 




+ h/FSV 
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+ X 


,} 


R531 • R534 R532 
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(99) 


For representative values of the parameters, equations (90) and (99) have 
been plotted in figs. 9 and 10. 


Case (ii) r'__/r = 0 and L'/L 4 0 ; (b) Uniform Combustion. 

S S 1 ' ~ r — ~ " " ■ ■ ■ ■ 1 ■ — 

The only difference between this case and that of adiabatic fluctuations 

j 

(treated under (a)) is that the wall temperature gradient is given by equation 
(74) and not by equation (96). Now we get 

(i + A + h 

' r-w t t / T— \ i ' 
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( 100 ) 


in place of equation (98). 

Using the convenient notation, 


R535 s (l )/ X 

N FSV ' ' 


we get 


t' = R535 — - 2 • R535 — 

W r= — 

P r 


( 101 ) 


Combining equations (95) and (101) to eliminate t' , we get 


w 


o _ R535 - (R532/R531 ) 

(R5 33/R53 1 ) + 2 • R535 


( 102 ) 
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Discussion on the Predicted Response Functions 

Several features of the present study invite a discussion on the derived 
response functions in order to bring into sharp focus their distinct differences 
from the response functions of previous theories. 

For our purposes here, the response functions of contemporary 
theories may be taken as adequately represented by 

(103) 


ft = 


~ - ( 1+A) + B 
1 X ] _ 


where A and B (not to be confused with the A and B of the nomenclature of 
the present report) are parameters related, to the activation energy for the 
surface pyrolysis reaction. For any given n for a propellant, different 
curves are obtained by the choice of different values of A and B, a pro- 
cedure that essentially implies having a flexible value of the activation energy 
for the pyrolysis reaction. In the present study, a single value of the acti- 
vation energy (28. 9 k cal/mole) has been used throughout. That is, the fa- 
miliar A-B parametric representation of the response functions is outside 
the scope of the present work. The propellant is allowed to determine its own 
pressure index n through the reaction rate and the vapor pressure criterion, 
as contrasted with other theories where the pressure index n is specified as 
an empirical number from external considerations. The proper zero frequency 
(\j -* 1) limit of ft (i. e. , ft^ _ ^ = n ) is automatically reached in those theories, 
for the simple reason that this limit is used as one of the requirements on the 
response function in the derivations. In the present study, no such require- 
ment is imposed, and hence some interesting observations can be made on the 
behavior of the response functions at very low frequencies. Two other features 
of the present calculations which were mentioned earlier are worth repetition. 
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The nature of the present theoretical analyses introduces a high frequency 
limit on the validity of the derived response functions. Since simplifications 
have been introduced by limiting our attention to moderate frequencies, the 
interpretation of the theoretical results above a non-dimensional frequency 
(0) of 10 is questionable. 

For the purposes of heat transfer calculations, the assumption of a 
homogeneous solid has been made all through the present study. That is, the 
assumption of the characteristic linear dimension (effective thermal depth, 
x/r) associated with the heat transfer process in the solid being far larger 
than the characteristic linear scale (oxidizer particle size, a) associated with 
solid heterogeneity is inherent in all of the present theoretical derivations. 
Since the thermal diffusivity may be taken as a constant during the small vari- 
ations in the wall temperature encountered, the assumption of homogeneous 
solid becomes questionable at high regression rates and hence at high pres- 
sures. This is an important consideration. Interestingly enough, results of 
unpublished experiments currently near completion have strongly indicated 
that the familiar stabilizing effect of high pressures on unstable burning pro- 
pellants may be due to the fact that the propellant is taken beyond the homoge- 
neous solid limit, so that the heat reservoir effect in the solid (i. e. , the 
C * AT "charge" - "discharge" effect) would lose its significance. Deprived 
of the mechanism introducing phase differences in heat -transfer related vari- 
ables, the propellant would burn in a non-oscillatory manner above a certain 
"threshold" pressure. 

Following the above general comments we now examine the specific 
cases presented in figs. 8, 8b, 9, 10, and 10b. 



- 66 - 


Figure 8: Deep Solid Reactions with No Surface Reactions; Adiabatic 
Fluctuations in the Gas Phase. 

For this particular case, the results have already been discussed in 
Section 5. 2 . The essential point to note is that the "flame" temperature is 
assumed to vary with pressure according to the adiabatic law, an assumption 
that is difficult to justify at very low frequencies and is definitely incorrect 
at zero frequency (as revealed by thermochemical equilibrium calculations). 
The present analysis relates the wall temperature, and hence the regression 

i ' 

rate through the pyrolysis law, to the flame temperature. Thus, increases 

in flame temperature that accompany a pressure increase lead to large in- 

! 

creases in the regression rate. However, if the flame temperature does not 
vary at low frequencies, neither does the wall temperature, and hence we ap- 
proach the steady-state pressure index n.i The fact that no peaks are revealed 
in the curves at any frequency indicates that the values of the parameters used 
(E/RT w , in particular) do not permit such a behavior. The mean pressure de- 
pendence of the response function is due to the vapor pressure equilibrium ef- 
fect. When wall temperature fluctuations are present, the fragment size and 
the pyrolysis rate are not affected in the same manner as during time- 
independent combustion. This is essentially because of the phase differences 
(between the pressure and the wall temperature fluctuations), which are in- 
operative during time -independent variations in mean pressure. For example, 
if the flame temperature fluctuations vanish, the regression rate variations at 
low frequencies are affected by pressure variations exactly as during steady 
combustion, and hence the response function tends to the limit of steady-state 
pressure index n for all values of the mean pressure. 
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Figure 8b: Deep Solid Reactions with No Surface Reactions; Uniform 
Combustion in the Gas Phase . 

i 

This case provides a particularly interesting example in the wake of 
the discussions on the adiabatic case (fig. 8). For the case of quasi- steady 
uniform combustion in the gas phase, the gas-phase temperature variations 
accompanying pressure fluctuations are almost the same as those during time- 
independent combustion. However, at higher frequencies, the phase lag in- 
troduced by the condensed phase heat reservoir effect does not affect the 

pressure sensitivity in the reaction rate term and the vapor pressure term in 

| 

the same way, so that the curves reveal a mean pressure dependence. 

Figure 9: Surface Reactions in a Melt Layer of Constant Thickness 
Augment the Deep Solid Reactions; Adiabatic Fluctuations in the Gas Phase. 

When the thickness of the melt layer on the surface does not vary 
during oscillatory burning, all of the oscillatory burning characteristics are 
controlled by the deep solid reaction behavior. Thus, we see that the results 
(i. e. , frequency dependence of the response function) are likely to be close to 
the case of no surface reactions, as indeed they are in fig. 8b. The rest of 
the discussion on fig. 9 would be very similar to the one on fig. 8 and is not 
pursued here. 

Figure 10; Surface Reactions in a Melt Layer of Fluctuating Thickness; 
Fluctuating Part of Regression Rate Associated with Deep Solid Reactions 
Vanish; Adiabatic Fluctuations in the Gas Phase . 

The propellant studied is the one for which the time -independent char- 
acteristics were presented in figures 5 and 6. In the analysis (Section 5. 3) it 
was assumed that the regression rate fluctuations resulted from variations in 
melt layer thickness only. That is, accompanying a pressure increase there 
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would always be present a decrease in the melt-layer thickness of a magnitude 

precisely requisite to balance out the regression rate variations (see eq. 83). 

; 

If such variations occurred during time -independent combustion, the propel- 
lant would exhibit a zero pressure index, i as shown in fig. 5. It would be a 
simple matter to consider an incomplete balance and thereby introduce non- 
zero values of the pressure index n. That is, a propellant with any n can be 
considered almost as easily. The zero-n case represents one interesting 

■ i 

limit in the general class of propellants. 

Turning our attention to the frequency dependence, we expect such 
propellants to be fairly unstable from physical considerations. The deep solid 

I 

I 

is homogeneous and gives rise to temperature fluctuations due to the heat 
reservoir effects. The surface reactions rates are hence subject to large 
fluctuations because of the strong Arrhenius term. It should be recalled that 
for the present case of constant-temperature melt layer, regression rate de- 
pends directly on the surface reaction rate and not On its square root. Visual- 
izing the phenomenon in another way, the temperature fluctuations originating 
from the deep solid effects drive the Arrhenius rate term at the surface. The 
pressure dependence of the position of the peak on the frequency scale is also 
easily understood. With the existence of a melt layer on the surface, we have 
a new length scale in the problem, the thickness of the melt layer (£). We 
thus form a natural frequency, r/ l. For a constant mean regression rate 
with the variation of mean pressure (i. e. , the zero-n propellant) the melt 
layer thickness decreases with increasing pressure (see fig. 5). Hence, the 
natural frequency increases with increase of mean pressure. If we associate 
some "resonance-like" phenomenon with the peak in the response function, it 
would seem physical to expect such a peak to occur close to this natural fre- 
quency. This explains the shift towards higher frequencies of the peak with 
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increases in pressure. 

Figure 10b: Surface Reactions in a Melt Layer of Fluctuating Thick- 
ness; Fluctuating Part of Regression Rate Associated with Deep Solid Reactions 

Vanish; Uniform Combustion in the Gas Phase. 

. 

For the case of a zero-n propellant considered here, the response 
function is not expected to differ appreciably from the one with adiabatic 
fluctuations in the gas phase. As can be seen, most of the discussions on the 
previous figure are independent of the nature of the temperature or pressure 
fluctuations in the gas phase. The broader peak observed in fig. 10b is 
probably due to much milder temperature fluctuations compared to the adia- 


batic case. 
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5.4 Further Considerations in Gas Phase Processes 

Several aspects of the problem provide a strong motivation for the 
study of gas phase in detail. It was seen in Section IV that a completely 
self-contained solution to the time-independent combustion problem has not 
been possible because of our lack of understanding of gas phase details. 

In the present section we observed that our progress toward a solution to 
the oscillatory combustion case was also hindered by the same difficulty; 
more specifically, by our inability to write a general expression for the 
flame temperature fluctuation as a function of the pressure fluctuations and 
frequency. In many of the experiments on oscillatory combustion, the 
pressure fluctuations are seen not to grow indefinitely but to reach a limit- 
ing value, indicating the presence of nonlinear effects. However, the pres- 
sure fluctuations themselves are almost purely sinusoidal. Thus, the 
origin of the nonlinearity is not obvious. It has been shown recently that 
under typical conditions the rocket chamber cannot sustain nonlinear pres- 
sure oscillations so that the nonlinearities are associated with either the 
condensed phase or the gas phase combustion zone. The condensed phase 
behavior was seen (Secs. 5. 2, 5. 3) to be controlled strongly by the gas 
phase "flame region" through the relation T^ = fn(p', id, . . . ). The above 
thoughts add to the fundamental interest intrinsic with the unconventional non 
premixed "flame" by itself. 
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From dimensional considerations we write 


T -T 

,, , b w , , ua XV > 

W = 1ST » “ 


w o 

Since ® ~ V , ua/$ may be considered as a Reynolds number, 
thus have 


R. We 


W = fn(R, — ) . 

a 

Except at very low values of R (low r and hence low p ), molecular dif- 
fusive effects are not likely to be important. Hence, at sufficiently high 
burning rates, 


W = fn( — ) 
a 

independent of gas velocity. This is reminiscent of the classical experi- 
ments of Hawthorne, et al. on combusting turbulent jets where the brush 
length was a function of feed velocity in the low flow-rate regime only. At 
high flow rates, the brush length was independent of feed velocity. The 
added complication in the case of propellant combustion is that the "feed" 
velocity (u) is dependent on wall temperature (through the pyrolysis law) 
so that the group x*/ a would be a constant only if W is constant. Intuitively 
one might expect for all propellants the existence of such an asymptotic 
limit at which both the wall temperature and the flame standoff distance be- 
come truly independent of regression rate. The numerical values of the 
limits, however, cannot be determined from dimensional considerations 
alone. A formal solution would require a quantitative knowledge of the mix- 
ing processes. We also recognize that equation (26) is no longer adequate 
to describe the flow, unless the property values, k , p , and c are aver- 
aged over macroscopic transport also. 
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During oscillatory combustion, the pressure waves in the far field 
(chamber) are acoustic. Thus, the temperature fluctuations are known as 
a function of the pressure fluctuations. This is valid up to the downstream 
end of the combustion zone, e (fig. 1). We need the amplitude of the tem- 
perature fluctuations at the upstream end of the combustion zone, b (fig. 1 ). 

It is suggested that an experiment be performed to gain insight into 
the problem. To facilitate measurements, the composite propellant sur- 
face would be simulated by a perforated porous plate with provision for in- 
jection of fuel gas through the fine mesh and oxidizer gas through the per- 
forations. The fraction of total surface area occupied by the perforations 
will simulate the oxidizer volume fraction in the propellant. The scaling 
rules for feed velocities are given by the foregoing dimensioned, considera- 
tions. The choice of gaseous fuel and oxidizer gases would control the 
flame temperature T^ . Measurements of wall temperature and the flame 
standoff distance at various pressures would complete the time -independent 
combustion studies. 

The same apparatus could then be used to study the case of oscilla- 
tory burning. The same details as before would be studied with isentropic 
pressure disturbances in the main chamber. The coupling between the 
fluctuations of feed rate and the temperature and pressure is difficult to 
simulate precisely. However, we are mainly interested in the behavior of 
the acoustic wave across the nonlaminar combustion zone. The experi- 
ments are expected to shed light on an important facet of the general prob- 
lem of composite propellant combustion. 
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VI. CONCLUDING REMARKS 
Most of the significant developments have already been dis- 
cussed in the previous sections. However, the following points call for 
special mention. 

/ (i) From dimensional considerations one recognizes a 

natural regression rate r= n/a for all composite propellants. Interesting 
characteristics are anticipated for propellants operating in the vicinity of 
this natural regression rate. At values of r < r , the propellant may be 
looked upon as homogeneous and hence deep solid reactions may be im- 
portant. At values of r > r , the thermal profile in the solid is too shal- 
low to lend credence to the homogeneous - solid approximation. Surface 
reactions are expected to dominate the regression rate behavior. The 
present study indicates (Section IV) that the pressure index n should 
change by approximately a factor of two as one passes through the r = r 
point. Many propellants have been experimentally found to exhibit such a 
pronounced "break" in the r versus P curve (T-17, CIT-2 . .. for ex- 
ample). If the oxidizer particle melts and mixes thoroughly with the 
readily melting binder, the particle size a loses its significance and 
hence there exists no natural regression rate. Such propellants are not 
expected to show a marked break in the r versus P curve. 

(ii) If, by some means, the degradation of the oxidizer is 
suppressed in the subsurface region, the pressure index n is expected 

to be unusually high for a composite propellant. This offers a simple 

25 

explanation for the high n obtained by Muzzy when he coated the oxi- 
dizer particles with a passive layer. Also, the surface reactions, with 
little deep solid reactions, could give a very pronounced peaky response 


G 
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function as shown through the present analysis (Section 5. 3 and fig. 10). His 

25 

experimental data provide a strong support for the theoretical developments 
here. 

(iii) One of the most important developments of the present work is 
the identification of the pressure index n as a composite quantity which re- 
flects the combined influences of several different physical phenomena. Thus, 
a zero-n propellant can exhibit strong instability, if these different influences 
are not affected in a self -compensating way during oscillatory burning (Sec- 
tion 5. 3 and fig. 10). This offers an explanation for the behavior of zero-n 
propellants whose tendency to oscillate has not been explained so far on the 
basis of existing theories which treat n as a simple entity completely charac- 
terizing the pressure sensitivity of a propellant both during steady and un- 
steady combustion. The identification of the different physical influences be- 
hind the composite entity n also shows how the response function can be de- 
pendent on mean chamber pressure, contrary to previous theoretical results, 
but in agreement with experimental data. 

(iv) It should be recognized that the heat of combustion (Q, Q^) is pres- 
sure dependent for a given propellant. As defined here, the heat of combustion 
represents the heat released in the gas phase by vapors leaving the propellant 
surface. The net heat release would depend on the initial endothermicity of the 
propellant vapors. The initial endothermicity depends on the number of un- 
broken bonds left in the fragments entering the vapor phase. The number of 
unbroken bonds is characterized by the fragment size vaporizing FSV. For 
vapors of a specific fragment size (which may be taken as the reference value), 

A 

the heat of combustion would have a fixed value (Q ). For vapors of mean initial 
fragment size larger than this reference value, the heat of combustion would 
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be smaller than Q, since the larger number of unbroken backbone bonds left 
in the fragments introduce a larger heat sink on the total heat release. The 
reference state is arbitrary and may be taken as unity FSV. Now, the heat 
of combustion of propellant vapors is, 

Q f = Q - D(1 - ~V ) ‘ (104) 

Recalling that the mean fragment size vaporizing (FSV) depends on the wall 
temperature and the chamber pressure through the vapor pressure equilibrium 
criterion, we see that the heat of combustion also depends on the wall 
temperature and the chamber pressure. During time -independent burning, 
such variations in are expected to slightly modify the results of the uni- 
form combustion model in this report. Preliminary calculations have shown 
that the modifications are no more than a few per cent in r and hardly any at 
all in . Likewise, during oscillatory burning, we expect fluctuations in 
the value of because of fluctuations in the value of FSV through equation 
(104). The response functions (employing the uniform combustion model) in 
the present report are expected to undergo slight modifications because of this 
variable value of Q^. The results in the report have used a constant value of 
. Again, it is expected that the modifications would be minor; nevertheless, 
the variable case is anticipated to remove the small discrepancy in the zero 
frequency limit of the response function in fig. 8b. A variable was con- 
sidered in a crude manner in Section 5. 2 (page 55). A more rigorous procedure 
is under consideration. 

(v) In the present report, a linear structure has been assumed for AP 
in order to relate the mean fragment size at any state (of degradation) to the 
number of unbroken backbone bonds in AP. This is a debatable assumption. 

A more accurate procedure would be to bring out this relationship through 
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the actual structure of partially degraded (but undocomposed) AP if the struc- 
ture is accurately known. The results of such a refined procedure are, how- 
ever, not expected to alter the present conclusions significantly, mainly be- 
cause the concepts of FSV and vapor pressure equilibrium at the wall would 
not be affected, and most of the principal conclusions of the present report 
result more from the concepts of FSV and vapor pressure equilibrium than 
from the actual value of FSV as related to the degradation rate. 

Salient Conclusions 

(i) The assumption of condensed-phase degradation of AP being the 
rate-limiting step in the overall pyrolysis of AP leads to theoretical results 
that match well with experimental hot-plate data. 

(ii) The hypothesis of pressure-dependent (condensed-phase) degrada- 
tion of AP in a surface layer controlling the overall regression rates of AP/ 
composite propellants leads to predictions in agreement with experimental 
data covering both single-crystal AP self- deflagration and composite propellant 
combustion. 

(iii) The extent of degradation at the vaporization step can be specified 
through the vapor pressure equilibrium criterion. 

(iv) The hypothesis of gas -phase combustion rate completely deter- 
mined by pressure-independent diffusive -mixing processes leads to analytical 
results in agreement with experimental trends. 

(v) The response function of a composite propellant can be theoretically 
derived through the method of inner and outer expansions including explicitly 
the non-linear (Arrhenius) degradation rate term in the condensed phase. 

(vi) The response functions so derived exhibit dependence on mean 
chamber pressure, strongly so with the model of adiatatic fluctuations in the 
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gas phase and weakly with the model of uniform combustion in the gas phase. 

(vii) A "zero-n" propellant cam exhibit fairly strong instability be- 
havior. 

(viii) A few unifying concepts enable us to coherently interpret a host 
of superficially diverse data. 

In conclusion, the author wishes to thank Professor Fred E. C. Culick 
for many helpful comments and encouragement, Mr. Warren L. Dowler for 
his interest in these studies and encouragement, and Mrs. Roberta I. Duffy 
for her excellent typing. 
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APPENDIX A. Details of Algebra 


(i) Equation (62) from Equation (61) 
We had 

r = (h+D^njj— + e 


w 


w h 


( 61 ) 


Using equation (54) to write £ in terms of £ , we write equation (61) as 

W 


X = (h+l)-fn 




h + 6° 


FSV 


+ h iL (\ 1 n , i 

FSV - V " X, / + 1 


1 J 
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Recognizing that 


FSV 

"S° + 1 r« / 


FSV- 1 

Lit + fsv - \ 

1 " X,/J 


we may approximate the logarithmic term and write 
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FSV 


Defining -X: = (h+1 ) <tn { j 
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+ 1 r ' (\ 
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+ FSV - V ' 
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■ — (] 
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, where FSV 

FSV 


(A3) 


the time -dependent components FSV and FSV' , equation (A3) is written, 
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Substituting the value of 


from equation (51) into equation (A4), 
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(62) 


(ii) Equation (63) from Equation (62) 

= (h+ 1 ) tn { FS y_ 1 } - j^gy 


During oscillatory combustion, 

FSV = FSV + FSV' . 


I 

Hence, 
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Approximating the logarithmic term, as usual, 
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From Appendix B it is seen that a good representation of FSV is 

-P ' 


FSV = a P" K exp(6T ) 
i w 


where a, (3, and 6 are constants. It is easily seen that 


i 

t-h nr* i I r r 

FSV' P' , w t a P' , w| w“ o K 

— - _ _Q + 0 = -p + O 
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-P — + t' 6X . (A6) 
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w wow 

Substitution of equation (A5) into equation (62) and utilization of equation (A6) 
yields equation (63). 


(iii) Solution of Equation (12) 

The following solution obtained through the method of matched 

asymptotic expansions is reproduced from reference 1. It would aid the 

reader to go over the descriptive discussion of the solution method on page 17. 

pp' + p = (h+p+T) A exp (-E/ftT) (12) 

T = 0:p = 0; t = 1 : p = p (13) 

w 


Outer solution 


Away from the wall region, the exponentially small re 

I 

(R. H. S. ) can be neglected from equation (12). 


14 + 41 = 0 
a7 dy 


T = 0 : p = 0 


action rate term 


( A7 ) 

(A8) 


i 
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p = -t is the solution to equations (A7) and (A8) to all orders. The solution 
cannot satisfy the hot boundary condition because equation (A7) is not valid 
near the hot boundary. 

Inner solution 

The large activation energy encountered in degradation reactions 
limits the temperature range of interest to very small values around the 
wall temperature. Hence, equation (12) may be approximated near the wall 
as 


pp'+p = (h+p+T)A exp {- 9 [l+X(l-T)l] . (A9) 

• cl 

Since A -is a very large number (it has B in it), it is convenient to define 

k - A exp( - 9 ) 
a 

and get 


pp'+p = (h+p+T )9r exp{-0 X(1~ T )} 

a 


(A10) 


Identifying the small parameter, e = 1/9 X , and defining the inner variables, 

q = ( 1 - t )/ e , £ = (p+t) , 

along with a new nondimens ional regression rate X = e-X , we write equa- 
tion (A10) as 

-(^-1+eTi)^' = (h+£ ) X exp(-rt) . (All) 

Here, a prime denotes differentiation with respect to r| . 

I 

At T = 1 (on the wall), r> = 0 


•n = 0 : £ =p+T = - 1 - — + 1 = _ . — 

w r w w FSV FSV 


Making the expansions, 


e = e D + e ^i + + e 3 ^3 + ••• • 

2 3 

X = X + gX . + e X- + s X- + . . . 
o 1 2 3 


(A12) 


and substituting these expansions into equation (All), we get to zero order 



/ 



Kaplun's matching theorem requires 

(inner limit of the outer solution) = (outer limit of the inner solution) 
or 

outer inner 
P t-» 1 ~ P rroo ’ 

which gives % - 0 as r| -* oo , yielding through equation (A15) the value of 
the constant as 

(h+ 1 ) tn h 

Substituting this value into equation (A15), we get 

(h+1 K,n(h+£ )-£ = -Xexp(-ri)+ (h+1 )hih . (A16) 

Using the boundary condition on the wall, namely, 

^ = • Wv ° n 11 = ° » 

we get 

X = (h+1 Kn{ psy.j} - psv » 
which was quoted earlier as equation (14). 


( A.l 7 ) 
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APPENDIX B. The Fragment Size Vaporizing FSV 

The importance of the statistical mean size (FSV) of fragments 
vaporizing from the wall cannot be overstated. Specification of FSV through 
a scientific procedure (which has often been taken as arbitrary in the older 
literature) was seen to hold the key to a host of problems in hybrid com- 
bustion. Throughout the present study we have had ample opportunities for 
recognizing the role of FSV in the deflagration of composites. We thus feel 
the need for a quick rule for evaluating the FSV as a function of wall tempera- 
ture and chamber pressure. 

After some study of API hydrocarbon vapor-pressure data, the fol- 
lowing rule was evolved to predict the fragment weight as a function of inter- 
face temperature (°K) and chamber pressure (atmospheres). 

W. = 32. 8 P" 0 ' 2615 • exp(3. 67 x 10" 3 T ) . (Bl) 

The equation has been evolved to slide rule accuracy and no minimum- 
error procedure has been employed. The equation has not been tested beyond 
the temperature range available in API charts (i. e. , T = 1200°F or 922°K). 
However, the extremely smooth data trend within this range appears to 
justify extrapolations as have been necessitated on occasions in the applica- 
tion of the above rule to composite propellant burning. 

For AP (molecular weight 117. 5) the above rule takes the form 

FSV = 0.2795 p- 0 * 2615 exp(3. 67 x 10' 3 T ) . (B2) 

w 

It seems more than probable that the weak power law for pressure 
effect is an approximation for a logarithmic law (which may be evolved 
through more accurate curve fit procedures). Such a logarithmic law 
would also seem more physical. In any case, the rule has been found to 
possess very good accuracy. 

0 ■> 
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Consider a linear chain of (n +1) elements (beads) having n 

oo oo 

bonds in the chain. We follow the degradation of this chain to smaller frag- 
ments. If n^ bonds are broken in the original chain, the number of frag- 
ments originating is (n^ +1). The mean size of the fragments, i. e. , the 
number of elements in a fragment- chain is 


FS = 


n + 1 
oo 

n, + 1 
b 


The number of bonds broken in the original chain is the original number n 

minus the total number of bonds n^ remaining now. Hence, 

n + 1 
__ GO 

FS = n - n + 1 • 

oo r 


oo 


If n^ ~ 10 we may neglect unity in comparison with n^ , and write 


n 
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CD 


n - n 
oo r 
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1 - 
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oo 


Denoting by N the fractional number of bonds, n^/n^ , we get 
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N = 1 


FS V 


(B3 ) 
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APPENDIX C . THE MELT LAYER 

To eliminate any arbitrariness associated with the useful concept of 
a melt layer on the surface of regressing AP, the following treatment is 
given as a first approximation. This treatment is an almost exact repro- 
duction from ref. 8. 

Although we normally use in it the wall values of temperature and 
fragment size, equation (14) is quite general and predicts the regression rate 
of any homogeneous material if the quantities, temperature and fragment 
size can be specified at a plane in the material. A parametric plot of equa- 
tion (14) would look like fig. A-l. The higher the mean molecular weight 
(fragment size), the higher will be the regression rate. This is because a 
higher fragment size implies a smaller number of bonds to be broken from 
the deep solid value, and it can be accomplished at a higher rate. The 
important point to note is that the regression rate is the same for all planes 
in the material. 

3 5 

A clear discussion in the Encyclopedia of Chemical Technology 
indicates a correlation of the following form between the molecular weight of 
a substance and its melting point: 

log M = A + B T 

o m 

35 

Typical values of A q and B are available for families of molecules. In the 
degrading AP, the "molecular" weight at any plane is simply the product of 
fragment size at that plane and the molecular weight of the fundamental 
NH^C-f-O^ (117. 5). Hence, the melting point vs. molecular weight correla- 
tion would be of the form indicated by the broken line in fig. A-l. 

At a specified regression rate r^ let us follow the changes in a plane 
as the plane moves from the low-temperature interior to the high -temperature 
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surface. At an intermediate station 1 , the temperature is and the 

fragment size is FS ^ - At the intermediate station 2 , the fragment size 

is FS^ and the temperature is • Now, as shown in fig. A-l, this also 
happens to be the melting point of a material of molecular size FS^ • The 
AP melts at this temperature, and the rest of the travel of this plane till 
the wall, i. e. , temperature , is through the melt layer. That is, the 
intercept 2-3 represents the melt layer. Now consider a higher re- 
gression rate r^ . We can see through the above arguments that the melt 
layer now represents a smaller intercept 2-3 at B than at A . At the 
limiting regression rate r^ , the melting point of the fragment size at the 
wall happens to exactly equal the wall temperature, so that the limiting con- 
dition of zero melt-layer thickness is reached. At a still higher regression 
rate r^ , the melt layer is imaginary because the melting point of the frag- 
ment size at the wall is higher than the wall temperature. 

The temperature intercept 2-3 may be related to the physical 
thickness through the temperature profile solution to the energy equation. 

Although the limiting case of vanishing melt layer is exactly valid, it 
should be carefully noted that the preceding arguments are not applicable in 
a rigorous manner to the postulated liquid layer invoked elsewhere in this 
report. This is because the melt layers have been taken as isothermal in 
those applications; that is, different from the temperature solution to the 
energy equation as used here in Appendix C. It is expected that the argu- 
ments would, however, be valid in an "average" manner, considering the 
very small thickness of the melt layers encountered. 





Figure 1(c). Details of the Postulated Model for Composite Propellant Combustion. 



(a) Subsurface Reactions Only (b) Surface Reactions in a Melt Layer 

Augmenting Subsurface Reactions. 



Region of Active Degradation Reactions 

Approximate Order -of-Magnitude of Thermal Wavelength during Oscillatory Combustion 
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COMBUSTION OF AP-BASED PROPELLANTS 



Linear pyrolysis of ammonium perchlorate. Rot-plate pyrolysis: 
A Andersen and Chaiken (Ref. 7); O Coates (Ref. 8); □ ONERA (Ref. 1). 
Diffusion-flame pyrolysis (all from this work) : X in 760 mm Hg of (Hj + N 2 ) 
See Note 1; X in 30-35 mm Hg of (Hi + Ni) see Note 1; + in 100-200 mm 
Hg of (CH« + N») see note 2. 

Notes: 

(1) 1-in. diameter disks of AP burning in counterflow of gas with infrared 
surface temperature measurement, seo Ref. 3, Fig. lb(c). 

(2) Half-inch-square section rods, burning in parallel flow of (CHi + Nj) ; 
Ref. 18. 


2 . Theoretical Predictions for AP and Comparisons with 
Hot Plate Data from Powling9. 








Theoretical Predictions for a Typical Composite Propellant 
Considering Subsurface Reactions Only 

Mass fractions: 0.75 AP, 0.25 CTPB (binder kinetics from ref. 24) 
AP particle size, a: 90 (a 

Interfacial reaction layer thickness, \|l : 5 yj 
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Figure 4(a). 





Linear Regression in -l Rate, r (cm/sec) 










Typical Composite Propellant. Surface Reactions in Melt Layer Augmenting Subsurface Reactions. 



Melt Layer Thickness, 
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